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Both  the  cash  flow  and  the  timing  of  a  cash  flow  are 
often  found  to  be  random  variables  in  economic  analysis. 

Very  few  analytical  tools  exist  to  handle  such  problems. 

This  research  develops  analytical  tools  to  aid  in  the  analysis 
of  such  problems.  A  total  of  sixteen  basic  models  are 
developed  for  impulsive,  single  cash  flows  and  series  cash 
flows  involving  capital  investment  problems  with  mutually 
exclusive  alternatives.  Each  model  is  developed  under  the 
assumption  of  deterministic,  constant  interest  rates.  In 
addition  to  the  development  of  models  where  cash  flow  and 
timing  of  cash  flow  are  independent,  dependent  models  are 
also  presented.  The  use  of  these  model  elements  by  an 
analyst  will  allow  the  utilization  of  available  data  in  a 
more  realistic  and  accurate  mode.  The  models  are  also  of 

significant  worth  in  verifying  simulation  output. 
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Other  areas  covered  in  this  research  are  a  spanning 
set  of  models  found  in  the  literature,  the  development  of 
a  taxonomic  structure  for  capital  investment  problems,  a 
proposed  network  logic  for  solving  economic  problems,  and 


applications  of  this  research. 


Accession  For 
ST  IS  ORAkI 
OTIC  TAB 
Unannounced 
Justification. 


By - 

Distribution/ _ _ _ 

Availability  Codes 
lAvail  and/or  j 
Dlst  I  Special  \  ipccto 


MODELING  AND  SOLUTION  TECHNIQUES  FOR 
ADVANCED  ECONOMIC  ANALYSIS 


by 

John  H.  Estes  IV 
has  been  approved 
November  1981 


APPROVED. 


ACCEPTED. 


ACKNOWLEDGMENTS 


I  am  indebted  to  my  chairman,  Dr.  William  C.  Moor,  for 
his  constant  encouragement  and  for  the  many  hours  of  critical 
discussion.  I  would  also  like  to  thank  Dr.  D.  A.  Rollier 
for  his  guidance  and  help  in  the  verification  of  several  of 
the  mathematical  elements  of  this  research.  I  am  grateful 
to  Dr.  D.  D.  Bedworth,  Dr.  W.  E.  Lewis,  and  Dr.  R.  L.  Smith 
for  their  help  in  the  preparation  of  this  dissertation  and 
for  their  continual  counsel,  encouragement,  and  assistance 
in  my  studies. 

A  very  special  thanks  is  offered  to  the  faculty  and 
staff  of  the  Industrial  and  Management  Systems  Engineering 
Department  for  their  support  during  the  time  that  I  have 
spent  at  Arizona  State  University. 

It  is  with  deep  gratitude  that  I  thank  the  United 
States  Air  Force  and  our  great  nation  for  the  environment, 
confidence,  and  trust  given  to  me  during  my  career. 

I  express  sincere  thanks  to  my  wife,  Leanne,  who  always 
said  "you  can"  when  I  had  my  doubts,  and  to  my  children, 

Tracy  and  Johnny,  for  their  many  days  of  patience  while 
their  father  worked  on  this  project.  My  parents  deserve  a 
special  thanks  for  their  confidence  during  my  educational 
and  career  endeavors. 


iv 


TABLE  OF  CONTENTS 


Page 

LIST  OF  TABLES . - .  X 

LIST  OF  FIGURES . xii 

Chapter 

1.  INTRODUCTION  .  1 

Statement  of  the  Problem  .  2 

Approach  to  the  Problem .  2 

Organization  of  the  Dissertation  .  .  4 

Background  for  Engineering  Economics  .....  5 

Summary  .........  .  9 

2.  BACKGROUND  FOR  ANALYSIS  MODELS  .  10 

MAPI  Method . 11 

Morris  Model  . . 13 

Bowman  and  Fetter  Models  .  14 

Alchian  Model  .  . . 15 

Reisman  Models  ..........  .  18 

Bernhard  Model  .  20 

Hillier  Models  .  21 

Canada  and  Wadsworth  Model  .  .....  24 

Summary  .......  .  25 


v 


Chapter 


3.  TAXONOMY  FOR  CAPITAL  EQUIPMENT  EXPENDITURES  .  .  27 

Cash  Flow  Classification  . . 27 

Impulsive  Cash  Flows  .  .......  31 

Series  Cash  Flows  . . 35 

Continuous  Cash  Flows  . . 4l 

Cash  Flow  Parameters . 43 

Time  Classification  .  ......  45 

Interest  Rate  Classification  .  48 

State-of-the-Art  ...............  53 

Short  Form  Classification . 57 

General  Model  .  ......  58 

Summary  ........  .  60 

4.  ANALYTIC  TOOLS . 61 

Nomenclature  .  . . 62 

Independent  Model  Elements  .  63 

Gid^is  ”  Independent  Model  Element  .  63 

CS(iTss  -  Independent  Model  Element . 72 

CfgTig  -  Independent  Model  Element  .  7^ 

CssTss  -  Independent  Model  Element  .  78 

GisTcd  *  Independent  Model  Element  .  80 

CssTCd  -  Independent  Model  Element  .....  81 

GisTcs  -  Independent  Model  Element . 83 

CssTcs  -  Independent  Model  Element  .....  85 

Summary  of  Independent  Model  Elements  ...  35 

vi 


Chapter  Page 

Extension  for  Finite  Cash  Flow  Streams  ...  87 

Dependent  Model  Elements  .  .  ....  93 

Gid^is  "  Dependent  Model  Elements . 95 

Cid^is  -  Dependent.  Ramp  Model  Element  .  .  95 
GidTis  -  Dependent,  Decay  Model  Element  .  .  96 


GidTis  -  Dependent, Growth  Model  Element  .  97 

GisTis  “  Dependent  Model  Elements  .  98 

GisTis  "  Dependent,  Ramp  Model  Element  .  .100 
Gis^is  -  Dependent,  Decay  Model  Element  .  101 

GisTis  "  Dependent, Growth  Model  Element  .  102 

GisTcd  “  Dependent  Model  Elements  .  104 

GisTcd  “  Dependent,  Ramp  Model  Element  .  .  104 
Gis^cd  “  Dependent,  Decay  Model  Element  .  105 

GisTcd  "  Dependent.Growth  Model  Element  .  105 

GisTcs  -  Dependent  Model  Elements  .  107 

GisTcs  -  Dependent,  Ramp  Model  Element  .  .  107 
GisTcs  "  Dependent,  Decay  Model  Element  .  107 
GisTcs  "  Dependent,  Growth  Model  Element  .  108 
Summary  of  Impulsive,  Dependent  Cash 


Flow  Elements . 108 

CSdTss  -  Dependent  Model  Elements  .  110 


CsdTSs  “  Dependent,  Ramp  Model  Element  .  .  110 
GsdTss  _  Dependent,  Decay  Model  Element  .  112 
Gsd^ss  “  Dependent,  Growth  Model  Element  .  11 3 


vii 


Chapter 


Page 

CggTSs  -  Dependent  Model  Elements  .  113 

CSSTSS  -  Dependent,  Ramp  Model  Element  .  .  115 

Cgs^ss  -  Dependent,  Decay  Model  Element  .  115 

Cgs^ss  -  Dependent,  Growth  Model  Element  .  116 

Css^cd  “  Dependent  Model  Elements  .  116 

Css^cd  ~  Dependent,  Ramp  Model  Element  .  .118 

GssTcd  ~  Dependent,  Decay  Model  Element  .  118 

^ss^cd  -  Dependent,  Growth  Model  Element  .  119 

CssTcs  -  Dependent  Model  Elements  .....  119 

Summary  of  Series  Dependent  Model  Elements  .  121 

Extension  of  CSSTSS  Dependent  Model 

Elements  for  Finite  Cash  Flow  Streams  .  .121 


Expansion  of  General  Model  .  121 

5.  DISCUSSION  OF  MODELS . 125 

Taxonomic  Structure  .  125 

Independent  Models  .  126 

Dependent  Models  ........  .  128 

Summary  . . 129 

6.  NETWORK  REPRESENTATION  OF  A  CAPITAL 

INVESTMENT  ALTERNATIVE  .  136 

Present  Worth  Development  .  136 

Summary  . . 142 


I 


viii 


Chapter  Page 

7-  APPLICATIONS  FOR  RESEARCH  MODELS  .  143 

Extensions  to  Traditional  Engineering 

Economic  Problems  .  ....  143 

Capitalized-Cost  Comparisons  .  ...  144 

Example  7»1  •  . . 144 

Replacement  Economy . .  .  .  150 

Example  7.2 . 150 

Lease  or  Buy . 153 

Example  7*3  154 

Validation  and  Verification  of 

Simulation  Models  .  158 

Development  of  a  Simulation  Language  .  159 

8.  SUMMARY ,  SUGGESTIONS  FOR  FURTHER 

RESEARCH,  AND  CONCLUSIONS . 161 

Summary  . . .  161 

Suggestions  for  Further  Research  .  162 

Non-constant  Interest  Rate  Models  .  I63 

Alternate  Decision  Criteria  .  163 

Introduction  of  Additional  Parameters  .  .  .  164 

Simulation  Model  . . 165 

Conclusion . . . 167 

BIBLIOGRAPHY  .  . . 168 

APPENDIX 


I.  SELECTED  MATHEMATICAL  PROOFS  AND  DERIVATIONS  .  .  176 


IX 


Table  Page 

4.7*  GisTcd  -Dependent  Model  Elements . 106 

4.8.  C^sTcs  Dependent  Model  Elements  .  .  .  109 

4*9*  csd^ss  Dependent  Model  Elements  .  114 

4.10.  CSSTSS  Dependent  Model  Elements  .  117 

4.11.  CssTccj  Dependent  Model  Elements . 120 

4.12.  CssTcci  Dependent,  Finite  n  Cash  Flow 

Model  Elements  . . 122 

5*1.  Taxonomic  Structure  and  Current 

State-of-the-Art  .......  .  127 

5*2.  Mapping  of  Current  Research  to  Table  $.1  (5a)  .  130 

5*3*  Mapping  of  Current  Research  to  Table  5.1  (5b)  .  131 

5-4.  Mapping  of  Current  Research  to  Table  5*1  (5c)  •  133 

5*5«  Mapping  of  Current  Research  to  Table  5.1  (5d)  .  134 

6.1.  Repair  Data  for  Lathe  1 . 140 


LIST  OF  FIGURES 


Figure  Page 

3.1.  Examples  of  Discrete  Cash  Flows . 29 

3.2.  Examples  of  Continuous  Cash  Flows . 30 

3«3»  Heaviside  Unit  Step  Function  . . 32 

3*4.  Cash  Flow  Taxonomic  Structure . 46 

3*5*  Timing  of  Cash  Flow  Taxonomic  Structure  ....  47 

3*6.  Constant  Interest  Rate . 49 

3.7.  Stepwise  Constant  Interest  Rate  .  49 

3.8.  Variable  Interest  Rate . 49 

3>9.  Taxonomic  Structure  for  Interest  Rate  .  51 

3.10.  Taxonomic  Structure  for  Capital 

Equipment  Expenditures  .  52 

3.11.  Taxonomic  Structure  for  Research  Models  ....  56 

4.1.  Ramp  Function  .  . . 93 

4.2.  Decay  Function . 9^ 

4.3.  Growth  Function . . . 94 

6.1.  Present  Worth  of  a  Single  Cash  Flow  . . 137 

6.2.  Present  Worth  of  a  Multiple  Cash  Flow . 137 

60*  Network  for  Repair  Evaluation . 141 

7.1.  Dam  Construction  .....  .  l46 

7.2.  Pole  Replacement  vs.  Stubbing . 152 

7.3.  Owned  vs*  Leased  Truck  155 

8.1.  Simulation  Flow  Model  .....  .  166 


xii 


Chapter  1 


INTRODUCTION 

This  research  is  concerned  with  modeling  and  solving 
specific  classes  of  long-term  capital  investment  problems 
in  engineering  economics  involving  mutually  exclusive 
alternatives  with  deterministic  and/or  stochastic  parameters. 
The  analysis  is  based  on  the  net  present  worth  for  each 
alternative. 

Investment  decisions  are  among  the  most  difficult  and 
important  decisions  with  which  management  must  cope. 

There  are  several  reasons  as  to  why  a  detailed  analysis  of 
possible  alternatives  has  such  a  significant  impact.  First, 
the  decisions  involve  large  capital  expenditures.  Second, 
the  decisions  normally  have  long  lasting  impact.  The 
concept  of  irreversibility  in  the  decision  can  cause  a  poor 
decision  to  literally  destroy  the  firm.  Many  capital 
acquisitions f costing  several  thousands  of  dollars  may  have 
only  scrap  value  if  they  fail  to  achieve  the  objectives  of 
the  organization.  In  addition,  the  decision  may  very  likely 
commit  the  firm  to  a  plan  that  will  last  several  years 
prior  to  reaching  a  point  where  a  new  course  of  action  can 
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be  determined.  Third,  and  the  true  justification  for  this 
research  effort,  the  alternative  actions  invariably  involve 
a  high  degree  of  uncertainty.  They  involve  estimates 
concerning  many  different  variables.  When  the  variance 
and  probabilistic  characteristics  of  these  variables  are 
not  considered  in  the  decision  process,  the  decisions  are 
based  upon  incomplete  information,  and  in  some  cases,  can 
be  incorrect. 


Statement  of  the  Problem 

The  requirement  to  which  this  research  is  addressed 
is  the  development  of  solution  techniques  for  solving 
complex  alternative  selection  problems  in  capital  investment. 
Such  techniques  must  be  able  to  handle  a  wide  variety  of 
complex  discrete  and  stochastic  parameters.  The  general 
orientation  of  the  research  is  to  define  a  taxonomy  for 
capital  investment  model  elements,  to  formulate  and  analyze 
each  type  of  model  element,  and  to  establish  solution 
procedures  to  solve  complex  problems  in  the  area  of  study. 

Approach  to  the  Problem 

The  approach  to  the  problem  of  alternative  selection 
among  capital  investments  will  involve  three  basic  phases. 

The  first  phase  will  be  to  develop  a  detailed  classification 
for  capital  investment  models.  The  second  phase  will  be 
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to  discuss  analytical  tools  for  solutions  to  capital  invest¬ 
ment  problems  and  the  third  phase  will  be  the  actual  solution 
procedures  to  follow  in  solving  a  complex  problem. 

During  phase  one  of  the  research,  representative 
models  found  in  the  field  of  economic  analysis  will  be 
reviewed.  The  review  is  not  meant  to  provide  an  exhaustive 
enumeration  of  models.  It  will  introduce  the  reader  to  the 
current  state-of-the-art  in  solving  capital  investment 
problems . 

Fleischer  and  Ward  (1977)  have  developed  a  multiple- 
classification  taxonomy  for  economic  analysis  which  they 
claim  to  embrace  all  models  in  the  field.  The  taxonomy  is 
based  upon  separate  descriptive  classifications  for  cash 
flow,  interest  rate  and  planning  horizon.  Using  this 
taxonomy,  a  more  detailed  classification  structure  will  be 
proposed.  This  classification  structure  will  act  as  a 
guide  to  illuminate  those  areas  where  additional  research 
is  required. 

Specific  models  will  be  developed  and  selected  for 
further  study  during  phase  one  of  the  research.  These 
models  will  be  limited  to  those  alternative  selection 
problems  involving  deterministic,  constant  interest  rates. 

The  cash  flow  and  time  parameters  will  be  allowed  to  take 
on  discrete  or  continuous  and  deterministic  or  stochastic 


characteristics.  In  addition,  the  area  of  dependence 
between  cash  flow  and  time  will  be  discussed  for  those 
models  using  ramp,  decay,  and  growth  functions. 
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The  second  phase  of  the  research  will  be  to  develop 
analytic  tools  to  solve  capital  investment  problems.  A 
network  illustration  of  a  problem  will  be  discussed  as 
an  aid  in  illustrating  the  interrelationships  between 
cash  flow,  time,  and  interest  rate.  The  solution  procedures 
will  assume  independent  parameters. 

Phase  three  of  the  research  will  be  to  solve  a  complex 
alternative  selection  problem  using  a  general  model  and  to 
show  how  network  analysis  can  be  used  as  a  tool  in  the 
solution  process. 

Organization  of  the  Dissertation 

The  remainder  of  Chapter  1  will  be  devoted  to  the 
background  of  engineering  economics.  Chapter  2  contains 
a  spanning  set  of  models  used  to  solve  capital  investment 
problems . 

Due  to  the  wide  variety  of  problems  that  can  be  viewed 
as  capital  investment  problems,  Chapter  3  will  present  a 
taxonomy  of  such  problems  and  a  general  model  development. 

In  addition,  Chapter  3  contains  a  discussion  of  the  problems 
addressed  in  this  research. 


s 


Analytical  tools  for  solving  capital  investment 
selection  problems  will  be  discussed  in  Chapter  4.  A 
given  alternative  will  be  evaluated  by  calculating  a 
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present  worth  and  variance  of  present  worth.  The  remainder 
of  Chapter  4  will  be  devoted  to  the  development  of  a 
detailed  model  for  solving  complex  capital  investment 
problems.  Chapter  5  will  be  a  discussion  of  results. 

The  development  and  justification  of  a  network 
representation  of  a  capital  investment  problem  will  be 
given  in  Chapter  6 . 

An  analysis  of  a  specific  group  of  alternatives  is 
presented  in  Chapter  7.  This  analysis  is  presented  on  the 
basis  of  expected  present  worth  and  variance  of  present 
worth  for  each  alternative. 

Conclusions  and  recommendations  are  summarized  in 
Chapter  8.  Selected  mathematical  proofs  and  derivations 
are  given  in  the  Appendices. 

Background  for  Engineering  Economics 

Representative  textbooks  by  Fabrycky  and  Thuesen  (1964), 
Morris  (I960),  Terborgh  (1949)»  and  Taylor  (1964)  give 
general  credit  to  Wellington,  Fish,  Goldman,  Grant,  and 
Dean  for  the  development  of  the  three  basic  methods  of 
engineering  economic  analysis  (Rate-of-retum,  Present  value, 
and  Annual  cost). 


■» 


Buck  (1975*81)  attributes  the  roots  of  engineering 
economios  to  the  economics  of  railroad  building  at  the  end 
of  the  nineteenth  century.  The  initial  work  in  this  area 
was  accomplished  by  Wellington  (1908),  a  civil  engineer, 
who  was  interested  in  the  structural  selection  aspects  of 
railroad  building.  He  assumed  that  railroads,  due  to  their 
monopolistic  nature,  could  never  be  destroyed.  Also,  he 
set  forth  the  principles  of  compounding  in  assuming  that 
railroads  were  built  for  future  expansion  as  well  as  present 
requirements.  Future  traffic  was  postulated  to  increase 
on  a  compound  interest  basis.  His  book  was  addressed  to 
engineers  designing  railroad  systems.  The  focal  point  of 
his  work  centered  on  the  selection  of  structures  based  upon 
capitalized  cost. 

Another  civil  engineer,  J.  C.  L.  Fish  (1923),  published 
one  of  the  first  books  on  engineering  economics.  He  believed 
that  the  central  concern  of  engineering  economy  should  be 
in  choice  of  investment,  rather  than  structure. 

Goldman  (1923)  introduced  a  comparative  value  concept 
using  compound  interest  calculations  similar  to  those  used 
in  textbooks  today  (Lesser,  1969)*  His  work  leaned  heavily 
on  the  use  of  compound  interest  calculations  to  compare 
different  alternatives. 

In  1930,  Grant  published  the  first  of  a  series  of  books 
on  engineering  economy  (Grant,  1930)*  For  ‘the  first  time. 


short-term  investments  were  studied  in  depth.  Grant  ( 1938 1 
362-363)  pointed  out  several  obstacles  to  engineering 
economy  including  problems  with  data  leading  to  "hunches" 
and  oversimplification.  Other  problems  of  bias  on  the  part 
of  the  estimator,  impractical  alternatives  and  legal 
obstacles  were  also  pointed  out.  Despite  these  warnings, 
engineering  economy  is  still  largely  an  exercise  in  financial 
arithmetic  with  serious  omission  of  the  real  problems  in 
analysis.  Grant's  comments  (1938*81)  concerning  safety  are 
noteworthy  in  that  they  also  apply  to  any  risk  in  the  decision 
process. 

There  is  no  correct  way  to  introduce  into  an 
economy  study  this  requirement  for  a  margin  of 
safety  before  undertaking  a  proposed  investment  in 
plant. 

His  suggestions  to  include  safety  in  analysis  are  byi  using 
conservative  data;  increasing  the  required  interest  rate; 
using  a  mandatory  payback  period;  or  recognizing  the  need 
of  safety  as  an  irreducible  element  in  the  final  choice 
between  alternatives.  These  suggestions  are  based  upon 
deterministic  changes  in  the  data.  As  stated  by  Arthur 
Lesser  (19691  111),  "Eugene  L.  Grant  can  truthfully  be 
called  the  father  of  engineering  economy". 

Rate-of-return  analysis  was  introduced  by  Joel  Dean  (195i)» 
A  great  deal  of  argument  has  been  noted  in  the  literature 
over  common  errors  and  solution  problems  with  this  method; 
however,  Taylor  (1964*125)  notes  that  it  is  preferred  by  top 
management . 
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The  early  authors  such  as  Wellington,  Fish,  Goldman, 
Grant,  and  Dean  used  data  from  mortality  tables  or  curves. 
Rather  than  using  the  probabalistic  characteristics  of 
these  data  sources,  deterministic  data  were  used  by  cal¬ 
culating  means  and  using  these  figures  as  if  certainty 
was  assumed.  The  models  also  required  that  the  data  be 
put  into  a  specified  format  by  whatever  "appropriate" 
methods  and  the  "appropriate"  algebraic  manipulations  be 
used  as  stated.  There  is  very  little  discussion  as  to  how 
to  handle  risk  and  uncertainty  other  than  that  these 
elements  must  be  considered  in  the  decision  process. 

Modem  textbooks  have  developed  some  of  the  finer 
details  of  analysis  by  injecting  risk  and  uncertainty , 
utility  theory,  and  mathematical  and  computer  models.  The 
list  of  references  for  this  information  is  voluminous j 
however,  specific  references  to  some  of  the  more  sophisti¬ 
cated  works  is  required.  DeGarmo  and  Canada  (1973),  Morris 
(1 977)*  Raiffa  (1970),  Bussey  (1978),  and  Reisman  and  Rao 
(1972)  are  some  of  the  works  which  provide  powerful  tools 
for  analysis. 

The  work  of  Reisman  and  Rao  (1972)  is  perhaps  the  most 
detailed  effort  to  allow  for  stochastic  extension  of  model 
elements.  Interest  rates  and  inflation  are  considered  in 
both  discrete  and  continuous  modes  and  can  take  on  random 
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patterns  (Fleischer,  1975»79)*  The  equations  for  specific 
cases  are  derived  and  stated  in  a  set  of  tables. 

Summary 

A  review  of  literature  points  to  the  concern  in  society 
for  economic  analysis.  The  models  developed  to  accomplish 
the  analysis  use  a  large  number  of  parameters.  Since  one 
of  the  objectives  of  this  research  is  to  develop  a  general 
model  for  economic  analysis,  the  next  chapter  will  review 
a  representative  set  of  models  and  model  parameters. 
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Chapter  2 


BACKGROUND  FOR  ANALYSIS  MODELS 

There  are  several  well-developed  models  used  in  capital 
investment  decisions.  Since  one  of  the  objectives  of  this 
research  is  to  develop  a  general  model  to  evaluate  capital 
investment  problems,  a  review  of  representative  models  is 
essential.  The  selected  models  are  not  an  exhaustive  set, 
but  research  indicates  that  they  are  a  spanning  set  of 
models  used  by  industry  and  represent  the  "current  state-of- 
the-art".  The  models  span  analysis  based  on  parameters 
which  are  deterministic  to  parameters  which  are  stochastic, 
from  models  which  have  only  a  few  parameters  of  interest 
(less  than  or  equal  to  five)  to  models  having  as  many  as 
thirty-one  parameters,  and  from  models  which  are  concerned 
with  single  item  analysis  to  models  which  consider  a  string 
of  replacements  for  the  equipment  under  analysis. 

A  review  of  these  models  will  point  out  the  parameters 
which  were  considered  for  inclusion  in  the  taxonomy  presented 
in  Chapter  3.  The  purpose  of  the  review  is  not  to  critique 
these  models,  but  should  orient  the  reader  to  the  fact  that 
many  gaps  are  left  in  the  analysis  of  a  problem  due  to 


inapplicability  of  a  known  model  or  the  requirement  to  make 
simplifying  assumptions  to  the  data  used  in  the  analysis. 

The  review  will  also  be  used  for  the  general  model  to  be 
introduced  at  the  end  of  Chapter  3. 

The  models  selected  for  review  are  the  Machinery  and 
Allied  Products  Institute  (MAPI)  Method,  Morris  Model, 

Bowman  and  Fetter  Models,  Alchian  Model,  Reisman  Models, 
Bernhard  Model,  Hillier  Models,  and  Canada/Wadsworth  Model. 
Each  model  is  presented  with  a  brief  discussion  of  the 
history  of  the  model  along  with  its  mathematical  formulation. 
The  applicability  of  each  model  and  assumptions  made  from 
the  parameters  are  also  covered. 

MAPI  Method 

The  MAPI  Method  was  first  published  in  1958  by  George 
Terborgh  (1958).  The  MAPI  procedures  produce  an  "urgency 
rating"  which  is  based  upon  an  after-tax  rate-of  return 
analysis  of  the  net  project  investment  (Terborgh,  1958*153)* 
The  basic  method  is  to* 

1.  Select  an  ownership  period  which  may  or  may  not 
be  equivalent  to  the  economic  life  of  the  proposed  asset. 

2.  The  existing  asset  takes  on  a  one  year  life. 

The  analysis  is  a  one-more-year  rate-of-retum  analysis  to 
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compare  cost  between  keeping  the  existing  asset  an  additional 
year  and  replacing  the  existing  asset  with  the  proposed 
asset  * 

3.  The  formulas  capitalize  future  sums  beyond  the 
first  year  at  8.2 There  are  four  formulas  used  which 
vary  only  by  the  pattern  of  equipment  "inferiority"  used. 

By  pattern  of  equipment  "inferiority",  the  MAPI  method  is 
evaluating  the  rate  of  accumulation  of  depreciation  and 
obsolescence  over  the  service  life.  This  pattern  is  assumed 
to  be  straight-line,  double-rate  declining-balance,  sum-of- 
digits,  or  expensing. 

4.  An  urgency  rating  is  calculated  as  a  percentage 
on  the  extra  investment  after  recovery  of  all  costs  of 
buying  the  proposed  asset  now  rather  than  keeping  the 
existing  asset  one  additional,  year. 

One  of  the  four  formulas  is» 


C  =  n(Qn-wn)(Q-l)2  -  (1-b)  P  [  (Qn-1)  -  n(Q-l)]  _  ( q-1  ) , 
nQn(Q-l)  -  (Qn-l) 

where 


P  *  w” 


1  -  w  +  py  +  (l-p)z 
1-b 
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Q  =  1  +  i  -bpy, 


n  =  service  life, 
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p  =  fixed  debt  ratio, 

y  =  fixed  interest  rate  on  debt  capital, 

2  =  fixed  rate-of-retum  on  equity  after  tax, 

b  =  fixed  income  tax  rate, 

i  =  fixed  capitalization  rate, 

w1*  =  salvage  as  a  decimal  of  original  cost, 

w  =  ratio  of  salvage  in  year  n+1  to  salvage  in  year  n, 

C  =  next-year-capital  consumption  expressed  as  a  ratio 
of  investment  P. 

This  formula  is  for  an  accumulated  depreciation 
following  the  sum-of-the-year ' s  digits  tax  procedures 
(Taylor,  1964:374),  The  actual  solution  to  the  formula  is 
graphically  solved  and  a  set  of  worksheets  is  required. 

The  model  is  totally  deterministic  and  assumes  discrete 
parameters. 

Morris  Model 

William  T.  Morris  (1964:220-221)  presents  a  generalized 

equipment  replacement  model: 

TG(N,  ,No,  « •  •  ,Nj.,  . . .  )  —  Iq  +  ^  ■ — 2J — 

K  0=1  (l+i)'3 

-  S°N1  +  +  y2  CNp 

(l+i)Nl  (l+i)Nl  j=l  (l+i)Ni+j 

-  SN1N2  +  lN2  +  £2  °N2J 
(l+i)Nl+N2  (l+i)Nl+N2  j=l  (l+i)Nl+N2+J 


* 


. . . , 
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where 


It  =  initial  investment  in  a  machine  purchased  at  the 
end  of  period  t, 

=  operation  and  maintenance  costs  for  the  3th  year 
J  of  a  machine  purchased  at  the  end  of  period  t, 

S. .  =  salvage  value  at  the  end  of  the  jth  year  of  use 
J  for  a  machine  purchased  at  the  end  of  period  t, 

N.  =  life  of  the  kth  machine  in  the  sequence  of 
K  replacements. 

The  model  is  deterministic  and  involves  only  discrete 
parameters  with  single  yearly  costs. 


Bowman  and  Fetter  Models 


Bowman  and  Fetter  (1957*3^7-375)  present  two  models. 
The  first  model  is  for  a  chain  of  machines  or  replacements 
using  continuous  cost  and  revenue  functions  with  continuous 
interests 


rT 

J  [R(t)  -  E(t)]e_itdt  -  B  +  S(T)e~iT 
r  T  . 

/  R(t)e~ltdt  -  B  +  S(t)e'lT 
JO 


(l+e-iT+e"2iT+. . . ) 


1-e-lT 


where 


V 

B 

T 

R(t) 


E(t) 

S(T) 


present  worth  of  the  series  of  investments, 

constant  initial  cost  for  each  investment, 

life  of  a  piece  of  equipment, 

revenue  function, 

expense  function, 

salvage  value  at  T, 

annual  interest  rate. 


i 
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The  second  model  uses  the  same  basic  parameters  in  a 

discrete  mode. 

T  R+-E+  S(T) 

V  =  f  ~ — 2-  -  B  +  — - 

_t=l  (l+i)t  (l+i)T. 


• (l+i)T 
(l+i)T-l 

where 


V  = 


t=l  (1+i)^ 


-  B  + 


S(T) 


(1+i) 


T 


1  + 


(l+i)T  (1+i) 


2T 


=  revenue  for  a  given  period, 

E^  =  expense  for  a  given  period. 

The  two  models  require  that  all  investments  be  identical 
with  respect  to  all  parameters. 


Alchian  Model 


The  Alchian  Model  is  composed  of  an  expression  for 
expenses  and  an  expression  for  revenues  (Alchian,  1952). 

E  [ax  +  B1(l-e_wl.t)]  e"rtdt  +  (C-C0e"dln)e”rn 

oo 

+  C(l-ke-dL)  Z  e"r^n+L^ 

3*1 


oo 


Ae-2(n+jL)  +  B(1.e-»t)uLi  e-r(n+jL+T)dT  +  „ 


+  E  f[ 

j=o  Jo  L 

R  =  fn  (P  e"x,t)e"rtdt 

Jo  1 

oo 

+  I  (Lf  P  +  Q(i-e-g(ntLJ))l  e-sT-r‘n+L3+T>dT, 

.1=0  Jo  L 


where 
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A  =  initial  annual  rate  of  operating  and  main¬ 
tenance  costs  of  new  item  now  available, 

A.  =  current  annual  rate  of  operating  and  main- 

1  tenance  costs  of  incumbent  item, 

A  +  B  =  limiting  annual  rate  of  operating  and  main¬ 
tenance  costs  (deterioration)  of  new  item, 

A.+B.  =  limiting  rate  for  operating  and  maintenance 
1  1  costs  of  incumbent  item, 

C  =  purchase  price  of  new  items, 

CQ  =  salvage  value  of  currently  used  item, 

d  =  coefficient  related  to  rate  of  decline  of 
salvage  value  of  new  item, 

dj^  =  coefficient  related  to  rate  of  decline  of 
salvage  value  of  current  item, 

g  =  coefficient  for  rate  at  which  P  +  Q  is  approached, 

k  =  percentage  of  C  remaining  as  salvage  value 

immediately  after  purchase  (1-k  =  loss  due  to 
acquisition  and  installation), 

L  =  projected  replacement  period  of  new  equipment, 

n  =  economic  life  of  currently  used  equipment, 

P  =  initial  annual  rate  of  revenue  of  new  item, 

=  current  annual  rate  of  revenue  of  current  item, 

P  +  Q  =  limit  approached  by  initial  rate  of  revenue 

(prices  of  services  and  technological  changes 
in  new  machines)  as  time  passes, 

r  =  rate  of  interest, 

s  =  coefficient  for  rate  at  which  annual  revenue  of 
a  new  machine  changes  with  age  of  a  machine, 

s,  =  coefficient  related  to  rate  of  change  of  annual 
revenue  of  current  item, 


*» 


A 


* 


i .. .  •». >  -> 
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u  =  coefficient  related  to  rate  at  which  B  changes 
because  of  new  technology  and  prices , 

w  =  coefficient  of  rate  at  which  operating  and 

maintenance  costs  of  any  new  item  approach  the 
limiting  rate  (because  of  deterioration), 

w.  =  coefficient  of  rate  at  which  A*  +  B1  is 
approached  (because  of  deterioration), 

z  =  coefficient  related  to  the  rate  at  which  the 
initial  operating  and  maintenance  costs,  A, 
fall  as  new  items  are  developed, 

C0(t)  =  turn-in  value  of  present  machine  t  years  from 
now , 

A- ( t )  =  rate  of  operating  and  maintenance  costs  of 
present  machine  in  year  t, 

D(T)  =  turn-in  value  of  successor  machines  at  age  T 
as  a  fraction  of  C, 

A(t,T)  =  rate  of  operating  and  maintenance  costs,  per 

year,  of  a  machine  purchased  at  time  t  for  the 
Tth  year  of  its  age, 

e-rt  -  present  value  of  one  dollar  t  years  hence  at 
rate  of  interest  r  (continuously  compounded), 

t  =  time  in  year-units, 

j  +  1  =  number  of  items  in  series  of  machines;  j  =  0 
for  the  current  or  first  machine, 

V«(t)  =  value  of  annual  rate  of  services  in  the  year  t 
for  an  existing  machine, 

V(t,T)  =  value  of  annual  services  in  the  Tth  year  of  life 
of  operation  of  a  machine  installed  n  +  L(j-l) 
years  from  now. 

The  goal  of  this  model  is  to  select  the  combination  of 
parameters  which  yields  the  maximum  difference  R  -  E.  This 
goal  is  to  optimize  the  present  worth  of  the  appropriate 
replacement  model  by  use  of  computer  generated  tables  or 
dynamic  programming. 
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Reisman  Models 

Arnold  Reisman  has  published  detailed  models  for 
capital  equipment  investments.  There  are  three  basic 
approaches  used.  The  first  approach  is  to  consider  the 
initial  purchase  plus  installation  cost  and  the  salvage 
value  for  each  replacement  as  a  discrete  event.  The  expense 
and  revenue  parameters  are  considered  to  be  continuous  and 
all  discounting  is  continuous.  The  model  by  Reisman  ( 1971 « 73) 
ist 


where 


T,  =  time  at  which  the  hth  replacement  item  is 
installed , 

B.  =  cost  of  the  jth  replacement  item  which  is  con- 
J  sidered  to  be  the  total  cost  of  purchase  plus 
installation, 

Sj  =  salvage  value  of  the  jth  replacement, 

Ei(t)  =  expense  function  for  the  jth  replacement  which 
J  is  assumed  to  be  a  continuous  function  of  time, 
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Ri(t)  =  revenue  function  for  the  jth  replacement  which 
J  is  assumed  to  be  a  continuous  function  of  time, 

n-1  =  total  replacements  for  the  study  period. 

The  next  model  developed  by  Reisman  (1971*76)  considers 

all  cash  flows  to  be  discrete,  but  the  compounding  remains 

continuous.  He  modifies  his  basic  model  further  by  noting 

that  the  jth  replacement  can  have  multiple  cash  flows. 

Both  revenue  and  expanse  cash  flows  are  allowed  to  have 

k  cash  flows  within  the  economic  life  of  the  jth  replacement. 

The  model  is» 


P  = 


■n-1  j 

I  Be‘r£Th 
lj=o  Bie  h=0 

r-n-1 


r-n-1 


y  -r  f  (Th+1) 

Lj.oSl<T*+1>*  h“° 


V  a-rITh(fE  -rIT(j+l),q] 

i_n  e  _4r„Ejpe  q-o  ) 


13=0  -  -  'p=o 

n-1 


L  j=0  ”  ’  >o"  V 

The  following  model  treats  all  cash  flows  and  dis¬ 
counting  as  discrete. 

p  -  1  Jj,  V1+i)  ^  &  -  Jbs3(I3+l,(1+i)  [*i£H 


r  k 


L  (1+i)  t"h?oh^x  X^U+i)  ^q50T(j+1),q^ 


rn-l 

t  (l+i) 

L  3=0 


[-1*4 

h=0 


Lp=0 
k 


Z*Jp<l+i>  t-0LT(3+1)^ 


p=0 


p 

l 

q=0 


1 


"•*  t*  >V»Vr.  -  *V\ 
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It  should  be  noted  that  k  cash  flows  are  allowed  for  both 
revenue  and  expense  functions  for  each  replacement. 

The  next  set  of  important  work  accomplished  by  Reisman 
was  published  jointly  with  Arza  K.  Rao  (1973).  The  mono¬ 
graph  is  by  far  too  extensive  to  review  in  this  paper.  The 
research  sets  up  seven  deterministic  models,  similar  in 
structure  to  the  before  mentioned  Reisman  models,  which  are 
used  frequently  in  engineering  economics :  compound  amount 
of  a  single  payment,  present  worth  of  a  single  payment, 
amount  of  an  annuity,  periodic  deposits  to  accumulate  a 
future  amount,  present  worth  of  an  annuity,  capital  recovery, 
and  present  worth  of  a  deferred  annuity.  Each  of  the  seven 
models  is  then  evaluated  for  discrete  and  continuous  com¬ 
pounding.  The  models  are  further  extended  by  allowing  the 
cash  flows,  rates  of  discounting  or  compounding,  and  rate 
of  inflation  to  be  random  variables  with  specific,  known 
distributions.  In  addition,  the  interest  rate  and  inflation 
rate  are  allowed  to  be  time  dependent.  The  timing  between 
cash  flows  are  assumed  to  be  deterministic  and  discrete. 

Bernhard  Model 

Richard  H.  Bernhard  (1962j19)  presents  a  basic  model 
for  a  "proposed  productive  investment”. 
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P  =  Q0  +  Ql  +  Q2  +  . . .  _ Qs _ 

(l+i^  (l+ijMi+ig)  (l+i1)(i+i2)...d+is) 

+  •••  + _ 9a _  » 

(l+i1)(l+i2)...(l+in) 

where 

P  =  present  worth, 

Qg  *  net  incremental  return  to  be  gained  at  the  end  of 
period  s.  (s  =  0,1,2, ... ,n) , 

ig  =  rate  of  interest  on  borrowing  or  lending  in  any 
quantity  during  period  s. 

The  basic  model  considers  all  cash  flows  as  lump  sum 
end-of-period.  The  deviation  from  traditional  models  for 
discrete  cash  flows  is  in  recognizing  that  the  interest  may 
vary  between  periods. 

Hillier  Models 

The  Hillier  Models  (1963*449)  consider  cash  flows  as 
random  variables.  The  specific  relationship  between  the 
random  variables  is  assumed  to  be  independent  or  perfectly 
correlated.  Fleisher  (1975*7?)  points  out  that  Hillier' s 
work  was  the  most  important  contribution  to  capital  budget¬ 
ing  in  the  1960's.  Although  specifically  designed  for 
capital  budgeting,  the  basic  models  handle  any  random 
variable  cash  flow  with  known  means  and  standard  deviations. 
The  means  and  variances  for  the  present  worth  of  the  cash 
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flow  are  then  handled  analytically.  The  present  worth  for 
an  investment  isi 
n  X, 


P  *  £  7 - ryj  • 

j=0  (l+i)J 

where 

=  random  variable  cash  flow  during  period  j  having 
'  mean^j  and  standard  deviation^, 

i  =  interest  rate. 

Assuming  that  the  cash  flows  are  all  independent,  the 
mean  for  the  present  worth  of  the  entire  cash  flow  ist 


n 


Mp  =  I 


j=0  (l+i)J  . 


The  corresponding  variance  ist 

0% 


n 


'2  *  z 


j=0  (l+i)2j  . 


These  observations  are  fairly  straightforward}  however, 

Hillier  next  develops  the  random  variable  X.  ast 

J 

xj  =  rJ  +  Zj‘l)  +  Z3(Z>  *  ■■■  Z3U>' 

where 

Yj  =  independent  cash  flow  in  period  j, 

Z,(k)  =  kth  distinct  cash  flow  which  is  perfectly 
J  correlated  with  the  corresponding  cash  flow 

in  other  periods. 
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Since  Hillier  assumes  that  the  random  variables  are 


close  to  normal,  the  present  worth  should  be  close  to 
normal  with 


n 


(l+i)1' 


(1+i)3 


« 


and 


d2  = 


n 

Jo 


Var  (Yj) 

(1+i 


m 

+  Jo 


n 

ij=o  L 


^Var  (zT^) 
(l+i)i 


There  are  two  deficiencies  in  the  models.  First,  the 
assumption  of  normality  may  be  questionable  in  many  cash 
flow  structures.  If  it  is  assumed  that  the  number  of  cash 
flows  is  large,  the  resulting  present  worth  will  still 
have  an  approximate  normal  distribution;  however,  if  the 
number  of  cash  flows  is  not  large,  then  only  limited 
statements  can  be  made  using  the  Tchebycheff  inequality. 

The  second  questionable  assumption  is  in  considering 
that  the  correlated  cash  flows  are  perfectly  correlated. 
The  basic  concept  of  cash  flows  being  completely  indepen¬ 
dent  or  perfectly  correlated  is  at  best  questionable 
(Hillier,  1963*449).  To  consider  cash  flows  which  fall 
somewhere  between  mutual  independence  and  complete  corre¬ 
lation  requires  a  covariance  matrix  with  proper  weighting 
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factors  for  each  term.  From  an  analytical  point  of  view, 
the  development  might  not  be  too  difficult;  however,  from 
a  practical  standpoint,  it  would  appear  to  be  difficult 
enough  to  obtain  good  values  for  means  and  standard 
deviations  and  not  very  realistic  to  obtain  covariance 
type  data  from  investment  analysts. 

Canada  and  Wadsworth  Model 

The  model  proposed  by  Canada  and  Wadsworth  (1968)  is 
PV($)  =  P  +  D(l-e"rT)/r  +  Se_rT, 
where 

PV($)  =  present  worth  of  the  cash  flow, 

P  =  initial  investment, 

S  =  salvage  value, 
r  =  nominal  continuous  interest. 

Since  D  is  a  constant  receipt  or  disbursement  for 
each  year,  the  basic  model  differs  from  classical  models 
only  in  using  a  uniform  series  present  worth  calculation. 
The  importance  of  the  Canada  and  Wadsworth  work  is  that 
they  extended  the  traditional  model  for  conditions  where 
two  variables  such  as  salvage  value  and  time  could  be 
dependent.  This  was  accomplished  by  allowing  the  variables 
S,  T,  D,  and  P  to  be  dependent.  If  any  two  of  the  variables, 
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say  S  and  T,  are  dependent  while  the  remaining  variables 
are  independent,  then  the  variance  of  PV($)  is  approximately 
equal  to 


+  [6pV($)  l2  V(S)  +  2|6pV($)\/6pV($)  \  Cov  (T,S), 
'S'  T  II  S  / 

where 

V  =  variance, 

(5  =  partial  derivative  operator. 

An  additional  terra  is  added  for  each  set  of  dependent 
variables. 


The  models  presented  in  this  chapter  illuminate  the 
large  number  of  interi elated  parameters  which  are  considered 
to  be  of  importance  in  economic  analysis.  Authors  have  a 
rather  significant  difference  of  opinion  as  to  which  set 
of  parameters  should  be  used.  Although  each  model  is 
correct  when  viewed  from  the  specific  assumptions  made  in 
the  modeling  environment,  departures  from  these  assumptions 
invalidate  the  model. 


The  next  chapter  will  be  devoted  to  examining  the 
parameters  used  in  economic  analysis  to  develop  a  taxonomic 
structure  to  classify  different  problems  or  models.  This 
effort  will  be  used  to  approach  the  modeling  problem  by 
examining  various  model  components  or  elements.  The  taxonomic 
structure  will  represent  an  orderly  method  of  presenting 
models  and  will  also  illuminate  those  areas  where  additional 
research  is  needed. 


Chapter  3 


TAXONOMY  FOR  CAPITAL  EQUIPMENT  EXPENDITURES 

The  three  basic  parameters  for  all  capital  equipment 
expenditures  are  cash  flow,  time,  and  interest  rate.  Each 
of  these  parameters  can  take  on  a  multitude  of  various 
characteristics  such  as  being  discrete  or  continuous, 
deterministic  or  stochastic,  independent  or  dependent,  etc. 

A  review  of  the  literature  indicates  that  the  majority  of 
these  parameters,  along  with  their  associated  characteristics, 
have  not  been  researched  to  the  point  of  model  development. 
This  chapter  will  be  devoted  to  the  design  of  a  taxonomic 
structure  for  capital  equipment  expenditure  problems. 

Two  major  benefits  are  derived  from  the  taxonomic 
structure.  First,  it  is  imperative  that  a  classification 
scheme  be  used  to  point  out  areas  where  research  is  limited 
or  non-existent.  Second,  the  taxonomic  structure  presented 
in  this  chapter  is  used  to  illustrate  various  models  in  a 
logical  or  orderly  format. 

Gash  Flow  Classification 

The  cash  flows  in  a  problem  are  not  limited  by  the 
classical  approaches  to  solving  engineering  economic 
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problems.  In  the  classical  approach,  one  considers  all 
cash  flows  to  be  of  the  same  class  (e.g.  deterministic, 
continuous,  etc.).  However,  in  the  real-world  problem, 
one  might  consider  some  costs  such  as  scheduled  maintenance 
to  be  deterministic,  lump  sum,  and  independent}  whereas, 
another  cost  such  as  electrical  consumption  might  be 
deterministic,  continuous,  and  dependent  upon  operating  time, 
while  a  revenue  of  salvage  value  might  be  stochastic,  lump 
sum,  and  time  dependent.  With  this  in  mind,  it  appears  that 
one  must  first  consider  a  cash  flo w  model  as  being  composed 
of  different  cash  flow  model  elements.  Cash  flow  model 
elements  can  be  dependent  upon  each  other,  as  well  as  de¬ 
pendent  upon  time.  When  one  is  working  with  dependent  cash 
flows,  a  knowledge  of  the  covariance  matrix  for  the  cash 
flows  is  required  to  calculate  the  variance  of  the  present 
worth.  If  one  can  assume  that  the  cash  flows  are  independent, 
the  computation  of  an  expected  present  worth  and  variance  of 
the  present  worth  is  significantly  simplified. 

It  is  also  important  to  distinguish  between  expenses 
and  revenues.  A  simple  sign  convention,  as  normally  used, 
takes  care  of  this  problem. 

For  simplicity,  one  might  consider  cash  flows  as  being 
either  continuous  or  discrete.  Figure  3*1  is  a  representation 
of  possible  discrete  cash  flows  and  Figure  3*2  represents 
possible  continuous  cash  flows.  A  continuous  cash  flow  as 
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used  in  this  research,  includes  piecewise  continuous  cash 
flow  functions.  By  piecewise  continuous,  it  is  meant  that 
within  a  time  interval  there  are  a  finite  number  of  dis¬ 
continuities  and  that  the  right  and  left  limits  for  each 
cash  flow  within  a  subinterval  exists.  For  example,  Figures 
3.2(b)  and  3.2(d)  represent  piecewise  continuous  cash  flows. 

When  discussing  discrete  cash  flows,  it  seems  relevant 
from  a  mathematical  standpoint  to  distinguish  an  impulsive 
cash  flow  from  a  series  cash  flow.  An  impulsive  cash  flow 
is  essentially  a  cash  flow  which  is  not  related  to  any  other 
cash  flow  (e.g.  independent).  Examples  of  impulsive  cash 
flows  might  be  the  initial  purchase  cost  or  the  salvage 
value  of  a  piece  of  equipment.  Obviously,  any  series  type 
cash  flow  composed  of  n  elements  could  be  handled  as  n 
impulse  functions.  However,  the  ease  of  mathematically 
handling  series  cash  flows  justifies  the  additional  effort 
in  developing  both  models. 

Impulsive  Cash  Flows 

Several  authors  such  as  Fleischer  and  Ward  (1977il4) 
and  Reisman  (1971)  argue  that  continuous  cash  flows  and 
continuous  discounting  should  be  used  for  analysis.  The 
concept  of  continuous  cash  flow  and  continuous  discounting 
arises  from  the  argument  that  expenses  and  earnings  are 
created  every  second,  minute,  hour  and  day  of  plant  oper¬ 
ations  . 


32 


Using  continuous  cash  flows  and  continuous  discounting 
requires  that  discrete  cash  flows  be  treated  as  impulse 
functions.  As  developed  by  Spiegel  ( 196? i 255-258) ,  the  unit 

impulse  function  (t-tQ)  is  defined  by  its  integral  property j 
r  oo 


/ F(t )  6  (t-a)dt  =  F(a), 


(3.1) 


-'-oo 

where  F(t)  is  any  function  that  is  continuous  at  t  =  a. 


If  F(t)  is  set  equal  to  1,  then 
co 


I  (5(t-a)dt  -  H(t-a) 


(3.2) 


-CO 


‘  1,  t  =»a 
■  0 >  t  *=  a  t 


where  H(t-a)  is  the  step  at  t  =  a  as  illustrated  in 
Figure  3 . 3  and  is  referred  to  as  the  Heaviside  unit  step 
function. 


H(t-a)  I 


r 

i 

i 

i 

j_ 

a 


■>  t 


Figure  3.3 

Heaviside  Unit  Step  Function 
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One  can  now  look  at  the  Laplace  transform  of  6(t-a)  as 
being  equal  to 


f°° 

I  e"rt  (5(t-a)dt  =  e"ra. 

Jo 


(3*3) 


By  Equation  (3.1),  one  can  formally  derive  an  expression 
using  a  cash  flow  (C)i 


.  00 


Ge"rt  (5(t-a)dt  *  Ce-ra. 


(3-4) 


Equation  (3*4)  is  the  familiar  form  for  the  present  worth 
of  a  single  future  cash  flow  occurring  at  time  *  a  and  with 

continuous  interest  rate  r. 

Table  3.1  shows  a  variety  of  selected  functions  and 

their  associated  present  worth  formulas.  These  functions 
have  been  used  by  Hill  and  Buck  ( 1974 t 121)  and  other  authors 
due  to  their  common  occurrence  in  cash  flow  problems.  The 
step  function  is  found  commonly  in  operating  costs*  the 
ramp  function  is  typically  used  with  maintenance  and  deter¬ 


ioration*  the  decay  function  can  be  associated  with  startup 
and  learning  costs,  and  the  growth  function  is  commonly 
found  in  wear-in  maintenance  costs.  These  formulas  are 
developed  under  the  assumption  of  continuous  compounding. 
They  may  be  derived  for  discrete  compounding  by  denoting  i 
to  be  the  effective  interest  rate  and  derivingi 


er-l 

.-rt 


(i+1)1 


(3.5) 


Table  3.1 


Selected  Single  Cash  Flows  and  Transforms 
for  Continuous  Compounding 


Type  F(t)  Present  Worth  (Laplace  Transform) 


Step 

C 

Ce-rt 

Ramp 

Ct 

Cte-rt 

Decay 

Ce*at 

ce-(a+r)t 

Growth 

C(l-e-at) 

Ce"r't(l-e_at) 

Legend* 

C  is  the  scale  factor 
t  is  the  time  interval 
r  is  the  continuous  interest  rate 
a  is  a  constant 
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A  simple  substitution  of  Equation  3*5  allows  the  develop¬ 
ment  of  the  present  worth  with  discrete  interest  as  pre¬ 
sented  in  Table  3.2. 

The  procedure  for  handling  impulse  functions  is 

rather  simple.  Select  the  appropriate  type  of  cash  flow 

which  is  representative  of  the  expense  or  revenue  being 

considered  and  multiply  by  e"r-t  or  1  The  justification 

(1+i)* 

for  this  effort  will  be  further  developed  under  continuous 
cash  flows. 

Series  Cash  Flows 

As  pointed  out  by  Taylor  ( 1964: 31 ) ,  disbursements 
for  some  equipment  increase  with  the  life  of  that  equipment. 
By  similar  logic,  other  disbursements  and/or  revenues  may 
increase  or  decrease  during  the  life  of  the  equipment. 
Typically,  the  cash  flows  are  placed  into  partitions  with 
signed  components  and  multiplied  by  the  appropriate  series 
factor.  These  series  factors  require  single  calculations 
for  each  time  series  rather  than  by  treating  each  point  in 
time  as  a  separate  calculation.  This  analytic  advantage 
increases  with  the  number  of  points  considered.  Traditional 
approaches  consider  cash  flows  which  are  uniform,  arith¬ 
metic,  or  geometric  series  starting  at  time  equal  to  zero. 


36 


Table  3.2 

Selected  Single  Cash  Flows  and  Transforms 
for  Discrete  Compounding 


Type 

F(t) 

Present  Worth 

Step 

C 

C 

(l+i)* 

Ramp 

Ct 

Ct 

(l+i)* 

Decay 

Ce~at 

Ce-at 

(1+i)^ 

Growth 

C(l-e-at) 

C(l-e“at) 

(l+i)* 
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Hill  and  Buck  (197^* 120-125)  have  expanded  the  traditional 
analytic  tools  by  using  Zeta  transforms  which  serve 
engineering  economic  analysis  of  discrete  time  series  much 
in  the  same  manner  as  the  Laplace  transforms  do  those  with 
continuous  cash  flow  functions. 

Consider  a  function  F(n)  which  describes  a  cash  flow 
at  time  n.  The  present  worth  of  this  cash  flow  at  an 
effective  interest  rate  i  ist 
F(n)(l+i)“n. 

This  assumes  that  i  has  been  adjusted  or  that  n  is  a 
multiple  of  the  interest  compounding  period.  Given  that 
F(n)  is  a  series  of  end -of -period  cash  flows ,  the  present 
worth  would  bet 

f  F(n)(l+i)"n, 
n=o 

where  F(N)  is  the  last  cash  flow.  Hill  and  Buck  (19?4i83) 
note  that  this  looks  like  the  Zeta  transform 
oo 

I  F(nT)(l+zT)"n 
n=0 

after  one  replaces  i  with  zf  considers  T«l, and  assumes  an 
infinite  series  of  cash  flows.  One  notes  that  Zeta  trans¬ 
forms  can  be  used  directly  to  obtain  present  worth 
calculations  for  any  series  having  a  Zeta  transform.  Since 
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Zeta  transforms  sure  not  commonly  found  in  most  handbooks, 
the  Z  transform  tables  can  be  used  by  noting  that 
Z  =  1+zT. 

Table  3.3  is  generated  by  using  Z  transforms  and  then 
converting  to  the  appropriate  Zeta  transformation  (letting 
i=z).  To  illustrate  this  procedure  the  Z  transform  for 
F(nT)  =  CnT  is 

GTZ 

(Z-l)2  . 

which  is  equal  to  the  Zeta  transform 
CT ( 1+iT )  ■  CT(l+iT) 

(1+iT-l)2  (iT)2  . 

Authors,  such  as  Taylor  ( 1964 » 26),  develop  series  cash 
flows  which  start  at  time  equal  to  one  rather  than  zero. 

To  convert  Table  3.3  to  formulas  for  cash  flows  starting 
at  time  equal  one,  the  cash  flow  function  need  be  only 
evaluated  at  n»0  and  subtracted  from  the  appropriate  Zeta 
transform  entry* 

For  discrete,  series  cash  flows  which  do  not  start  at 
time  equal  to  zero  and  do  not  continue  infinitely,  adjust¬ 
ments  have  been  derived  by  Hill  and  Buck  (1974i123)* 

Table  3.4  represents  series  cash  flows  which  are  translated 
forward  b  time  units,  start  at  time  h  (where  hib),  and 
stop  at  time  k-1  (where  k»h). 
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Table  3*4 

Present  Value  of  Shifted/Translated  Series 
Cash  Plows  with  Discrete  Interest  i 


Type  of 

Series  F(nT)  Zeta  Transform 


Step  C 

Ramp  Cn 

Decay  Ce-an 


C 

i 


(l+i)1_h  - 


(1+i)1 


-h 


(l+i)1-* 

(l+i)1“k 


(h-b)(l+i) 


1-h 


C 

i 


(k-b)(l+i)1’k 


~^ha 


-  e 


-ka  "| 


C(l+i)1~b- 
l+i-e~a  L(l+i)h  (l+i)k 


Growth 


C(l-e_an)  C(l+i)*-b 
l+i-e"a 


l+i-e-a~ie~h^- 

i(l+i)h 

-  H-i-e~a4-ie~ka 
i(l+i)k 


Source i 

Thomas  W.  Hill  and  James  R.  Buck.  "Zeta  Transforms, 
Present  Value,  and  Economic  Analysis."  AIIE  Transactions, 
6,  No.  2  (1974),  123. 


The  logical  extension  of  the  material  in  Table  3.4  is 
to  develop  the  same  discrete  series  cash  flows  under 
conditions  of  continuous  compounding.  By  observing  the 
relationship  between  i  and  r  in  Equation  3*5»  the  results 
displayed  in  Table  3.5  can  be  easily  derived. 

Continuous  Cash  Flows 

The  derivations  for  continuous  cash  flows  can  be 
theoretically  derived  with  both  continuous  and  discrete 
compounding  of  interest.  However,  when  discrete  interest 
is  used,  the  model  is  really  a  subset  of  discrete  cash 
flow  with  discrete  interest  since  the  integration  of  the 
cash  flow  function  can  be  treated  as  a  single  lump  sum 
element.  For  this  reason,  discrete  compounding  is  part 
of  the  taxonomy  only  for  the  need  of  completeness.  In 
application,  one  would  simply  integrate  the  cash  flow 
function  over  time  and  multiply  by 

1 

(l+i)t. 

A  continuous  cash  flow  starting  at  time  h  and  con¬ 
tinuing  to  time  k  is  a  function  of  time,  C(t).  This 
function  is  then  multiplied  by  e_rt  to  discount  the  cash 
flow  to  the  present  and  integrated  over  the  time  period 
to  develop  an  equivalent  present  value  asi 
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Table  3.5 

Present  Value  of  Shifted/Translated  Series 
Cash  Flows  with  Continuous  Interest  r 


Type  of 
Series 

Step 

Ramp 


Decay 


Growth 


F(nT) 

C 

Cn 


Ce"*11 

Cd-e"311) 


Zeta  Transform 


e  -1 
C 


er( 1-h)  _  er(l-k) 

er(  1-h)  -  er(l-k) 


(er-l)2  L 


+  _c_ 

er-l 


er-l 


(h-b)(erU-4)) 
(k-b)er^1'k^ 


Cer(l-b)  r e_h(a-r)  _  e-k(a-r) 


er-e~a 


Cer(l-b) 


er-e”^ 


er-e~a- ( er-l ) e~4a 


(er-l)e^ 

-  -er-e‘a+(er-1)e~ka 

(er~i)erk 
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Buck  (1975i81)  and  other  authors  give  credit  to 
Grubbstrom  (1967)  for  recognizing  Equation  3*6  as  a  Laplace 
transform.  This  is  a  significant  realization  since  a 
table  of  Laplace  transforms  found  in  many  mathematical 
handbooks  such  as  Selby  ( 19 75 * 506-515)  can  be  used  to 
convert  many  continuous  cash  flow  functions  to  a  present 
value  equivalent.  Table  3.6  is  an  extension  of  Table  3.5 
for  continuous  cash  flow  functions  with  continuous  com¬ 
pounding. 

The  use  of  Laplace,  2,  and  Zeta  transforms  in  analysis 
have  two  important  properties  when  available  for  use  in 
modeling  engineering  economic  problems.  One  of  the  two 
properties  is  the  linearity  under  addition  and  multipli¬ 
cation  by  a  constant  which  allows  one  to  take  transforms 
of  component  cash  flows,  scale  each,  and  then  sum  them. 
Another  useful  property  is  the  ability  to  shift  the  time 
pattern  forward.  A  number  of  transforms  have  been  derived 
and  are  available  in  several  references.  V’hen  these 
transforms  do  not  exist,  they  must  be  either  derived  or 
more  involved  calculations  must  be  made. 

Cash  Flow  Parameters 

In  addition  to  classiiying  each  flow  as  impulsive, 
series,  or  continuous  one  must  also  look  at  the  parameters 


44 


Table  3.6 

Present  Value  of  Shifted/Translated  Continuous 
Cash  Flow  with  Continuous  Interest  r 


Type  of 

Cash  Flow  C(t)  Present  Value  (Laplace)  Transform 


Step 

Ramp 


Ct 


C(e-rh-e"rk) 


C(e-rh-e-rk) 


(h-b)e~r^  -  (k~b)e"rk 


Decay 

Ce~at 

Ce"br  j 

a+j 

e  ' 

Growth 

C(l-e~a^) 

£  (,-»*  - 

ce-rb  r 

a+r 

4-(a+r)h  _  -(a+r)k 
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associated  with  the  actual  cash  flow.  In  the  model  elements 
presented  prior  to  this  point,  deterministic  cash  flows  have 
been  assumed.  This  assumption  is  usually  made  in  traditional 
analysis;  however,  the  actual  cash  flows  can  also  be  random 
processes.  The  assumption  that  cash  flows  can  be  random 
variables  has  been  incorporated  into  the  cash  flow  taxonomic 
structure  illustrated  in  Figure  3.4. 

When  a  specific  cash  flow  is  studied,  the  timing  of 
that  cash  flow  must  also  be  studied.  The  timing  of  a  cash 
flow  is  the  next  subject  of  consideration. 

Time  Classification 

Once  a  cash  flow  has  occurred,  the  duration  or  time 
of  occurrence  is  known.  However,  prior  to  occurrence,  the 
timing  of  the  cash  flow  may  be  either  deterministic  or 
stochastic.  Once  again,  discrete  timing  is  divided  into 
impulsive  timing  for  a  single  time  element  or  series  timing 
for  a  number  of  time  periods.  Figure  3*5  represents  the 
time  classification. 

After  the  cash  flow  and  the  timing  of  that  cash  flow 
have  been  determined,  the  decision  must  be  made  for  an 
appropriate  discount  rate.  This  is  the  third  element  to 
be  classified. 
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Figure  3.4 

Cash  Flow  Taxonomic  Structure 
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Figure  3.5 

Timing  of  Cash  Flow  Taxonomic  Structure 
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Interest  Rate  Classification 

Three  basic  interest  rate  classes  are  considered  by 
Fleisher  and  Ward  (1977i21).  These  are  constant,  stepwise 
constant,  and  variable  constant  interest  rates.  Three 
examples  of  possible  interest  rates  are  illustrated  in 
Figures  3*6,  3*7*  and  3*8* 

The  interest  rate  is  defined  over  the  entire  time 
interval  of  a  specific  cash  flow.  The  actual  discounting 
takes  place  at  the  end  of  a  unit  of  time  or  continuously 
throughout  the  period.  In  addition,  the  interest  rate  can 
be  either  deterministic  or  stochastic.  In  the  deterministic 
constant  case,  the  interest  rate  (i)  is  fixed  over  the 
entire  discounting  period.  In  the  stochastic  constant 
case,  the  interest  rate  (i)  is  a  random  variable  that, 
once  determined,  is  constant  over  the  entire  discounting 
period.  In  the  deterministic  stepwise  constant  case,  the 
rate  is  constant  between  each  of  a  finite  number  of  dis¬ 
continuities.  The  end  points  may  not  agree  with  the 
discounting  time  period  which  will  require  an  adjustment  of 
the  discount  rate.  The  stochastic  stepwise  constant  case 
has  a  constant  mean  for  each  time  interval  and  each  random 
variable  ilf  ig,  ...  in  may  be  independent  or  correlated 

with  one  another.  In  the  deterministic  variable  case, 
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This  is  a  continuous  time  process.  A  similar  expression  is 
used  for  the  stochastic  variable  case.  Figure  3.9  illustrates 
the  taxonomic  structure  for  the  interest  rate. 

Inflation  has  been  considered  by  several  authors  as  a 
factor  which  must  be  considered  within  an  analysis.  The  work 
of  Reisman  and  Rao  (1973)  develops  both  the  justification 
and  procedures  for  handling  inflation.  For  the  purpose  of  the 
taxonomic  structure,  it  is  assumed  that  the  discount  rate  must 
be  appropriately  modified  to  reflect  inflationary  factors 
and  does  not  have  inflation  as  a  separate  classification. 

The  taxonomies  in  Figures  3.4,  3.5,  and  3.9  must  be 
combined  to  visually  depict  the  taxonomy  for  capital 
equipment  expenditures.  This  has  been  accomplished  in 
Figure  3. 10.  A  particular  calculation  for  present  worth 
would  be  formulated  by  taking  three  paths  through  the 
network.  As  an  example,  a  traditional  engineering  economic 
problem  is* 

Cash  Flow  -•  expense,  impulsive,  deterministic, 

Time  -  impulsive,  deterministic. 


Interest  Rate  -  constant,  deterministic. 


5- 


Figure  3*9 

Taxonomic  Structure  for  Interest  Rate 
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Note  that  a  classification  distinction  has  not  been 
made  for  continuous  versus  discrete  interest  since  it  has 
been  shown  that  this  involves  a  simple  mathematical 
substitution.  In  addition,  the  subject  of  dependence  among 
the  elements  of  cash  flow,  time,,  and/or  interest  rate  has 
not  been  visually  depicted.  Rather  than  redrawing  Figure 
3.10,  note  that  a  heavy  dashed  line  could  be  placed  between 
any  of  the  three  basic  parameters  to  signify  dependence. 

The  path  through  the  impulsive/series  classification  must 
agree  for  both  cash  flow  and  time. 

State- of- the- Art 

The  current  state-of-the-art  for  economic  analysis 
problems  is  reviewed  by  Fleischer  and  Ward  (197?t24,  exhibit 
3).  Their  work  is  reproduced  in  Table  3.7.  An  extensive 
review  of  the  literature  since  1977  does  not  indicate  any 
additional  entries.  As  can  be  noted  from  Table  3.7.  a  great 
deal  of  research  is  still  open  in  the  field  of  engineering 
economics.  Expected  values  for  some  decision  criterion 
such  as  present  worth  (Young  and  Contreras,  1975)  have  been 
developed  for  some  of  the  missing  entries  in  Table  3.7; 
however,  the  corresponding  variances  have  not  been  developed. 
In  addition,  only  the  works  of  Hillier  (1963,  1966)  develop 
dependent  models. 
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Table  3.7 

Taxonomic  Structure  and  Current 
State-of-the-Art 


Cash  Flow 

Inters:  Rate 

Amount 

Timing 

Constant 

Stepwise 

Constant 

Variable 

How 

Cert 

When 

Cert 

Det 

Sto 

Det 

Sto 

Det 

Sto 

Impulsive 

Det 

Disc 

Det 

Cla 

a 

Cla 

1 

Cla 

Sto 

Cont 

Det 

Cla 

Cla 

Cla 

Sto 

2 

Sto 

Disc 

Det 

EBB 

a 

mm 

Sto 

Cont 

Det 

wm 

mm 

M 

Sto 

B 

■  i 

■  i 

■  I 

Continuous 

Det 

Disc 

Det 

Cla 

ll 

Cla 

m 

Sto 

Cont 

Det 

Cla 

Cla 

Cla 

Sto 

Sto 

Disc 

Det 

2 

Sto 

Cont 

Det 

2 

Sto 

1- Reisman  and  Rao  (1973) 

2- Ward  (1975) 

3-  Hillier  (19 69) 

4- Hillier  (1963) 


Cert-Certainty 

Det-Deterministic 

Sto-Stochastic 

Cla-Classical 

Cont-Continuous 

Disc-Discrete 


Source 1 


G.  A.  Fleischer  and  T.  L.  Ward.  "Classification  of 
Compound  Interest  Models  in  Econon ic  Analysis."  Engineering 
Economist,  (Fall,  1977 )» 24. 


Any  one  of  the  missing  entries  in  Table  3*7  is  worthy 
of  research.  The  wide  use  of  minimum  attractive  rate  of 
return  (MARK)  in  economic  analysis  indicates  that  the  area 
for  initial  research  should  be  restricted  to  the  missing 
entries  in  Table  3.7  under  constant,  deterministic  interest. 
Grant,  Ireson,  and  Leavenworth  (1976s97-100»189-195)  pointed 
out  the  value  in  using  MARR  in  that  the  purpose  of  an 
economic  study  is  to  determine  if  an  investment  should  be 
made.  They  suggested  that  this  be  done  on  the  basis  of 
determining  whether  the  investment  can  be  recovered  with 
at  least  a  stipulated  MARR.  For  these  reasons,  the  models 
selected  for  research  in  this  paper  were  restricted  to  the 
missing  entries  in  the  impulsive  cash  flow  section  of 
Table  3.7  under  constant,  deterministic  interest.  Normally, 
an  impulsive  cash  flow  is  considered  to  be  a  single  payment 
or  expense  3uch  as  salvage  value  or  initial  purchase  cost. 
This  research  expands  impulsive  cash  flow  research  to  in¬ 
clude  series  (multiple  impulsive)  cash  flows. 

In  view  of  the  developed  taxonomy,  the  models  for 
Table  3.7  are  illustrated  in  Figure  3*11  by  talcing  one 
branch  through  cash  flow  and  one  through  time.  The  path 
through  the  impulsive/series  level  must  agree  for  both 
cash  flow  and  time.  The  dashed  line  indicates  dependence 
between  time  and  cash  flow.  The  lack  of  a  dashed  line 
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would  indicate  independence.  If  we  assume  that  revenue 
and  expense  can  be  treated  with  an  appropriate  sign  con¬ 
vention,  the  specific  models  can  be  enumerated  with  a 
short  form  classification. 

Short  Form  Classification 

The  concepts  for  the  classification  scheme  have  been 
addressed  in  prior  sections.  A  short  form  classification 
scheme  is  developed  in  this  section  to  enumerate  those 
models  studied  in  the  research.  The  proposed  form  is« 

CXYTZY, 

where 

C  =  Cash  Flow 

T  =  Time 

- i  impulsive  cash  flow 

X  = 

. s  series  cash  flow 
-d  deterministic 

Y  = 

-s  stochastic 
-i  impulsive  time 

Z  -  s  series  time 

-c  continuous  time 

A  bar  placed  over  C  and  T  indicates  dependence  between 
cash  flow  and  time.  Notice  that  interest  rate  is  not 
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included  since  each  model  assumes  deterministic,  constant 
interest.  Using  this  classification,  the  specific  models 
to  be  researched  are  presented  in  Table  3.8  for  independent 
cases.  Since  these  models  are  repeated  for  dependent 
cases,  a  total  of  1 6  models  are  developed.  It  is  important 
to  note  that  each  model  is  a  cash  flow  entity  where  a 
model  for  the  actual  present  worth  of  a  particular  alter¬ 
native  could  be  composed  of  a  large  number  of  independent 
models  for  each  cash  flow  entity.  Each  of  these  models 
are  developed  in  the  next  chapter. 

General  Model 

The  general  model  for  a  present  worth  analysis  is 
stated  asi 

E(PV)  =  E(P)  +  E(E)  -  E(R)  -  E(S) , 

where 

E  =  Expected  value  operator 
P  =  Purchase  price 
E  =  Expenses 
R  =  Revenues 
S  =  Salvage  value. 

Any  one  of  the  elements  (P,  E,  R,  S)  may  or  may  not 
contain  multiple  cash  flows.  This  general  model  is 
expanded  in  the  next  chapter  to  offer  more  detail  to  the 
analyst . 


Table  3.8 

Independent  Parameter  Research  Problems 


Model 

Cash  Flow 

Time 

GidTis 

impulsive 

deterministic 

impulsive 

stochastic 

Gsd^ss 

series 

deterministic 

series 

stochastic 

CisTis 

impulsive 

stochastic 

impulsive 

stochastic 

Gss^ss 

series 

stochastic 

series 

stochastic 

GisTcd 

impulsive 

stochastic 

continuous 

deterministic 

GSSTCd 

series 

stochastic 

continuous 

deterministic 

°isTCS 

impulsive 

stochastic 

continuous 

stochastic 

GSSTCS 

series 

stochastic 

continuous 

stochastic 
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Summary 

The  taxonomic  structure  for  capital  equipment  expen¬ 
ditures  illustrated  in  Figure  3. 10  represents  an  orderly 
and  logical  format  to  present  models  in  the  following 
chapter.  Several  characteristics  of  the  parameters  have 
been  left  out.  A  more  detailed  structure  might  include  the 
actual  functions  being  used  for  each  parameter  and  con¬ 
siderations  for  inflation.  As  pointed  out  by  Fleischer  and 
Ward  (1977«28),  "A  truly  comprehensive  classification 
scheme  would  be  so  extraordinarily  complex  that  it  would 
be  virtually  without  utility." 

Chapter  4  presents  models  for  the  selected  taxons  in 
Table  3.8  and  concludes  with  a  more  detailed  model  for 
analysis  of  capital  equipment  expenditure  problems. 


Chapter  4 


ANALYTIC  TOOLS 

The  models  proposed  for  study  during  this  research 
assume  lump  sum  cash  flow  elements  which  are  either  impulsive 
(single  cash  flow)  or  series.  Analytic  tools  for  indepen¬ 
dent  and  dependent  models  are  presented  in  this  chapter 
with  a  concluding  section  on  a  general  model  development 
for  a  capital  investment  alternative  selection  problem. 

A  complete  model  for  a  specific  investment  alternative 
may  contain  multiple  impulsive  or  series  cash  flows, 
however,  each  cash  flow  element  is  assumed  to  be  independent 
of  all  other  cash  flow  elements.  This  assumption  is 
required  if  one  is  to  consider  time  as  a  random  variable. 

The  interest  rate  is  assumed  to  be  independent  of  both 
cash  flow  and  time.  In  addition,  the  interest  rate  is 
assumed  to  be  constant  over  the  time  horizon  and  determin¬ 
istic  in  nature  for  all  model  elements.  Although  models 
are  developed  where  cash  flow  is  dependent  on  time,  the 
nature  of  this  dependency  is  restricted  to  ramp,  decay, 
and  growth  functions. 

Expected  value  of  present  worth  and  variance  of 
present  worth  are  developed  for  all  of  the  independent 


model  elements,  and  expected  value  of  present  worth  and 
variance  of  present  worth  are  developed  for  all  impulsive, 
dependent  model  elements.  The  variance  for  the  CsdTsg 
dependent,  ramp  model  element  is  also  derived.  The  complexity 
of  the  expression  degrades  the  use  of  it  as  an  analytical 
tool.  For  this  reason,  the  remaining  series,  dependent 
models  do  not  have  a  variance  expression}  however,  the 
expected  present  worth  is  derived  for  all  of  the  series, 
dependent  model  elements. 

Due  to  the  large  number  of  discrete  and  continuous 
distributions,  a  complete  listing  of  analytical  tools  for 
all  distributions  would  be  too  extensive.  The  approach 
used  in  this  research  is  to  develop  general  formulas  and 
then  show  applications  v/ith  commonly  used  distributions. 

Nomenclature 

The  short  form  classification  system  developed  in 

Chapter  3  is  used  to  introduce  each  model.  The  symbol  T 

is  used  for  time,  and  t  is  used  for  a  specific  numerical 

value  of  time.  Similarly,  C  is  used  for  a  cash  flow 

designation,  and  c  as  a  specific  value.  The  symbols  jj.  and 
2 

O  are  used  for  mean  and  variance  respectively  and  are 


subscripted  with  c  or  t  for  time  or  cash  flow.  The  symbols 
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E  and  V  are  used  for  the  expected  value  operator  and 
variance  operator  respectively.  Present  worth  is  denoted 
by  PW.  Due  to  the  extensive  use  of  moment  generating 
functions,  M-j.  and  Mc  are  used  for  the  moment  generating 
functions  for  the  time  and  cash  flow  distributions. 

The  independent  models  are  presented  by  first  deriving 
the  impulsive  model  elements  and  then  extending 
that  model  for  series  model  elements.  This  format  was 
selected  as  a  logical  development  for  the  independent 
models.  Since  variance  expressions  are  not  developed  for 
the  series,  dependent  model  elements,  all  of  the  impulsive, 
dependent  model  elements  are  presented  prior  to  the  series, 
dependent  model  elements# 

The  last  section  in  this  chapter  contains  a  more 
detailed  general  model  which  extends  the  general  model 
presented  in  Chapter  3  to  include  the  model  elements 
presented  in  this  chapter. 

Independent  Model  Elements 

Cash  flow  elements  in  this  section  assume  all  cash 
flows  are  independent  of  time.  As  is  normally  done,  the 
present  is  considered  to  be  at  time  equal  to  zero. 

cidTis  “  Independent  Model  Element 

The  first  model  discussed  is  when  the  cash  flow 
is  impulsive  and  deterministic,  such  as  with  an 
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overhaul,  and  time  is  impulsive  and  stochastic.  If  one 
designates  C  as  some  future  cash  flow  (impulsive  and 
deterministic )  then  the  present  worth  of  C  is 
PW  =  C-e"rr, 


and  the  expected  present  worth,  due  to  independence  of 
cash  flow  and  time,  is 

E(PW)  =  C  •  E(e-ri?).  (4.1) 


Since  T  is  a  random  variable  with  a  given  probability  law, 
the  solution  to  Equation  4.1  is  easily  found  by 

oo 

C-E(e-rf?)  =  C  •  Y.  e”rtP(t)  (4.2) 

t=0 


where 

t  =  (0,1, 2, 3.  ...), 

P(t)  =  probability  of  T  taking  on  a  specific  value  t. 
Since 


I  e"rtP(t) 
t=0 


is  the  moment  generating  function  for  P(t),  Equation  4.2 
can  be  written  as 


E(PW)  =  CMt(-r),  (4.3) 

where  M^(-r)  designates  the  moment  generating  function  of 
P(t).  The  derivation  of  the  expected  present  worth  for 
this  model  was  derived  by  Young  and  Contreras  (1975)* 
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Table  4.1  represents  some  common  discrete  distributions 
and  their  respective  moment  generating  functions.  To  find 
the  expected  present  worth  for  a  C^T^g  model  element,  one 
merely  multiplies  the  cash  flow,  C,  by  the  appropriate 
moment  generating  function  for  the  time  distribution. 

Before  proceeding  to  the  development  for  the  variance 
of  the  present  worth  for  this  model,  an  important  property 
for  present  worth  analysis  is  requiredj 

E(e-«):e*rE(?)  (4.4) 

The  vast  majority  of  engineering  economic  analysis  assumes 
that  Equation  4.4  holds  as  a  strict  equality,  as  an 
approximate  equality,  or  simply  ignores  the  effects  of 
random  variables.  Obviously,  in  those  cases  where  the 
data  is  completely  deterministic,  Equation  4.4  does  hold 
as  a  strict  equality. 

The  proof  of  Equation  4.4  is  derived  by  a  Taylor 
series  for  e-"1^  expanded  about  the  point  T  and  is 

found  in  Appendix  I  (1).  Ignoring  the  subject  of  variance, 
Equation  4.4  points  out  where  a  present  worth  analysis  can 
be  in  error.  When  one  of  the  following  conditions  exist, 
the  present  worth  analysis  of  a  group  of  mutually  exclusive 
alternatives  could  lead  to  an  incorrect  selection. 
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Table  4.1 

Discrete  Probability  Laws  and  Their 


Moment  Generating  Functions 

Distribution  Probability  Law  Mt(-r) 

Binomial  P(t)  =  (t)ptqn“t.  t=0,l,...,n  (pe"r+q)n 

Geometric  P(t)  =  pqt-1#  t=l,2,3,...  pe~r 

0,  otherwise  l-qe-r 

Poisson  p(t)  =  t=0. 1.2.3,...  e  l(e_r-l) 

tl 


0,  otherwise 


6? 


1.  One  or  more  of  the  alternatives  involve 
random  variables.  Since  the  present  worth  is  underestimated 
for  an  alternative  involving  time  as  a  random  variable,  an 
incorrect  selection  is  possible. 

2.  A  budget  limit  for  the  selection  of  alter¬ 
natives  is  required.  Ignoring  the  idea  that  time  is  a 
random  variable  may  place  the  least  cost  item  under  the 
budget  limit  when  it  is  in  fact  over  the  budget  limit. 

The  significance  of  this  error  is  important  and  must 
be  evaluated.  As  an  example,  let  us  assume  that  the  time 
between  major  breakdowns  for  a  capital  investment  follows 
the  Poisson  distribution  with  mean  =  A.=  1  year.  Further¬ 
more,  let  us  assume  that  the  repair  cost  is  independent  of 
the  time  before  a  breakdown  and  is  equal  to  $10,000.  The 
nominal  interest  rate  is  10%,  The  expected  present  worth 
by  Equation  4.3  isi 

E(PW)  =  C-Mt(-r) 

=  $10,000  e^"’10"1^  $9,092.25. 

Using  E(PW)  =  C  «  e~r^,  leads  to  a  solution  of  $9048.37. 

This  is  a  .483??  error.  Obviously,  this  is  a  very  small 
error  and  should  not  significantly  effect  the  selection 
process  if  this  alternative  was  among  other  similar  capital 
investment  alternatives.  Now  let  us  assume  that  the 
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nominal  interest  is  allowed  to  increase.  The  percent 
error  can  be  obtained  fromi 

CMt(-r>  -  Ce-r^t,100> 


GMt(-r) 

For  the  Poisson  distribution,  the  error  is 


,  A.e“r  .  e  A(l-r) 
.Xe-r 


100. 


Table  4.2  illustrates  the  error  as  the  nominal  interest 
rate  increases.  The  question  of  what  is  a  significant 
error  is  rather  academic j  however,  the  error  does  increase 
rapidly  as  one  uses  nominal  interest  rates  over  10$.  The 
current  use  of  interest  rates  well  in  excess  of  10#  would 
support  the  requirement  that  random  variables  should  be 
taken  into  consideration  for  a  present  worth  analysis. 

The  variance  for  the  present  worth  of  the  CidTig 
independent  model  can  be  derived  from 
V(PW)  =  E(PW2)  -  [e(PW)]  2 


=  C2 


EU"1*^2  -  [EU-3^]  2j  . 

The  derivation  is  given  in  Appendix  I  (2)  and  is  shown  to 
reduce  to« 


V(PW) 


M. 


;t(-2r)  "  [Mt(-r>] 


Table  4.2 

Error  for  C^^ig  Independent  Model 
(Time  Poisson  Distributed 

With  A. "  1  year) 
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In  the  example  problem,  the  variance  for  the  present 
worth  would  bet 

V(PW)  «  C2[llt(-2r)  -  [Mt(-r)]  2 

=  ($10.000)2[e‘e'2(-10>-1>  -  [,<»-10-D]  2 
•  $752,043.9^* 

The  standard  deviation  ist 

(S  m  $867.20. 

Without  the  use  of  moment  generating  functions,  Poisson 
tables  (PARZEN,  1960t  444)  must  be  used.  A  series  of  17 
calculations  are  required  to  obtain  the  expected  present 
worth  and  a  series  of  9  calculations  are  required  to  obtain 
the  variance. 

As  stated  earlier,  one  must  allow  for  more  than  one 
cash  flow  (sill  independent)  for  a  given  alternative.  For 
alternatives  which  involve  model  elements  exclusively 

the  expected  present  worth  for  the  kth  alternative  is» 

where 

nk  =  number  of  impulsive,  deterministic  cash  flows 
for  the  k’th  alternative, 

Cjjj.  **  the  impulsive,  deterministic  cash  flow  for 
J  the  k^h  alternative, 

Mt  (-r)  ss  the  moment  generating  function  for  the  time 

parameter  associated  with  the  jth  cash  flow  of 
the  k'th  alternative. 
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The  selection  model  (based  on  mean  present  worth)  is: 
nk 

“jf  1  “tjk  ‘-r)- 

Alderfer  and  Bierman  (1970041)  suggest  the  logical 
argument  that  selection  among  alternatives  with  tied  mean 
present  worths  should  be  based  on  which  of  the  tied  alter¬ 
natives  has  a  smaller  variance  and  is  used  for  the  tie¬ 
breaker  in  this  research.  Variance  is  also  discussed  in 
Chapter  7  for  making  probability  statements  concerning  the 
present  worth  of  an  alternative. 

The  variance  for  the  present  worth  of  the  k *4  alter¬ 
native  is: 


^k 


V(PWk) 


■  I 

3=1 


'3k 


LMtJk 


(-2r)  -  [Mtjk(-r)] 


In  the  event  that  the  moment  generating  function  for 
a  distribution  does  not  exist,  calculations  using  Equation 
4.2  must  be  used  to  find  the  expected  present  worth.  The 
variance  of  the  present  worth  must  be  found  by: 


[  oo 

oo  \ 

V(PW)  =  c2 

1  e-2rt  F(t) 

lt=o 

-(I  e-rtP(t))2 

t=0  '  / 
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However,  this  would  be  a  rather  unusual  calculation  in 
that  the  moment  generating  functions  for  most  distributions 
do  exist  and  are  found  in  a  multitude  of  handbooks  and 
textbooks.  The  appropriate  interrelationships  between  the 
Z,  Zeta,  Laplace  transforms,  and  moment  generating  functions 
demonstrated  in  Chapter  3  allow  for  a  large  number  of 
distributions  to  be  handled  without  lengthy  calculations. 

The  remainder  of  this  research  assumes  that  the  moment 
generating  function,  or  an  appropriate  transform,  does 
exist . 

Alternative  selection  models  for  the  remainder  of  the 
model  elements  are  discussed  in  the  final  section  of  this 
chapter.  The  introduction  of  the  alternative  selection 
model  for  C^dT^s  model  elements  was  presented  at  this  point 
to  orient  the  reader  with  respect  to  a  decision  choice 
methodology  of  minimum  expected  present  worth  among 
alternatives . 

The  next  model  developed  is  an  extension  of  the  C^T^g 
model  for  cash  flow  elements  which  are  impulsive  and  series 
in  nature. 

CgdTss  -  Independent  Model  Element 

If  one  now  considers  a  series  type  cash  flow  which  Is 
deterministic  in  nature  and  has  an  impulsive,  stochastic 


.  .  >*  .■ 
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time  of  occurrence  for  each  cash  flow,  the  present  worth  is 


calculated  byi 

PW  =  Cle-rtl  +  c2e-r(tl+t2)  +  c3e-r(tl+t2+t3) 

+  ...  +  Cne^^l+^+V'-^n)  +  ...  . 

If  the  Ci's  are  different  in  magnitude,  then  one  must  solve 
for  the  present  worth  by  considering  each  cash  flow  as  a 
C’id^is  model  element.  However,  if  the  cj/s  are  a  constant 
function (denoted  as  C),  then  the  expected  present  worth  can 
be  shown  to  be  equal  to» 

E(PW)  =  C  Mt(-r) 

l-Mt(-r)  ,  (4.5) 

and  the  corresponding  variance  can  be  shown  to  be  equal  toi 


v(pw)  =  C2[Mt(-2r)  -  (Mt(-r)]2] 
[l-Mt(-2r)]  [l-Mt(-r)p 


(4.6) 


The  derivations  for  Equations  4.5  and  4.6  are  quite 
lengthy  and  are  presented  in  Appendix  I  (3). 

As  an  example  of  the  use  of  Equations  4.5  and  4.6, 
let  one  assume  that  a  preventive  maintenance  program  costs 
$100.  each  time  maintenance  is  done.  The  time  between 
maintenance  is  assumed  to  be  a  function  of  running  time, 
which  in  turn  is  a  function  of  demand,  the  time  between 
maintenance  requirements  is  assumed  to  be  Poisson  distributed 
with  a  mean  =A.=  1  month.  The  nominal  interest  rate  is  24# 
(2#  per  month). 
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By  Equation  4.5#  the  expected  present  worth  isi 

e(pw)  -  loo  (e*e  ~1}  =  $5,000.33, 

(e(e-°2_1) 

V(PW)  a  (100)2  e(e~,0,Zf-1)  _  (e^e*",02~1^)  2] 

[l  -  e<e-°4-l>]  [l  -  e<e“02-^]2 
=  $255,008.66  , 
and  the  standard  deviation  is 

(S  =  504.98. 

Using  the  traditional  formulas  (Taylor,  1964i26)  with 
the  timing  assumed  to  be  deterministic  (one  month  intervals) 
result  in  a  present  worth  of  $5,000.  The  insignificant 
error  of  33  cents  is  expected  with  the  low  monthly  interest 
rates;  however,  having  the  variance  gives  one  more  infor¬ 
mation  for  a  decision.  This  is  illustrated  in  Chapter  7* 

CisT,a  -  Independent  Model  Element 

When  both  cash  flow  and  time  are  impulsive  and 
stochastic,  one  can  calculate  the  expected  present  worth 
(assuming  independence)  from: 

PW  =  C  e"rt  , 

E(PW)  =  E(C)  *  E(e~rt).  (4.7) 

Since  C  is  not  limited  to  being  discrete,  one  must 
develop  formulas  for  continuous  cash  flow.  Table  4.1 
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suffices  for  discrete  cash  flows  by  merely  substituting 
Mc(-r)  for  Mt(-r).  Table  4.3  represents  the  moment  gener¬ 
ating  functions  for  a  few  common  continuous  distributions. 

The  mean  and  variance,  as  developed  by  Whitehouse 
(1973«24),  for  continuous  distributions  arei 
E(C )  «/ie  -  dMe(-r) 

d(-r)  -r=0 

E(C2)  =  d2  Mc(-r) 

d(-r)2  -r=0 

Cfc2  =  E  (C-^J2  =  E(C2)  -  [e(C)]  2 
=  d%  (-r)  _  /dM  (-r)\2  I 


d(-r): 


d(-r) 


|-r=0 


In  terms  of  moment  generating  functions,  Equation  4.7 


can  be  expressed  asi 


E ( PW )  =  dMc(-r) 


d(-r )  -r=0 


Mt(-r) 


(4.8) 


Equation  4.8  can  be  written  simply  as 

E(PW)  =/ic  Mt(-r)  .  (4.9) 

The  variance  for  the  present  worth  follows  directly 
from  the  classical  variance  of  the  product  of  two  inde¬ 
pendent  random  variables  (Duncan,  1965»94)i 

V(PW)  -yLtc20e_rt2  +  K-rt)2tfc2+  <tfe-rt)2tf02,  (U-10> 


! 
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Table  4.3 

Continuous  Probability  Laws  and  Their 
Moment  Generating  Functions 


Distribution 

Probability  Law 

He(-r) 

Exponential 

f(t)  =  ae“at,  t-0 

a 

0,  t«0 

a+r 

Gamma 

f(t)  =  a  (at)k“*eat, 

P(bj  0, 

t-0 

t-=  0 

l±f 

Normal 

f(t)  =  1  e-^l±=H)2 

(5W  1 

-  oo«=  t  •*=  1 

9 

DO 

e-rj^r2CTc2 

Uniform 

f("t)  =  1  ,  a  *=  t  <b 

-rb  -ra 

e  -  e 

b-a 

-r(b-a) 

0,  otherwise 
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(1  =  Mean  cost  or  revenue 


fJLe- rt  =  Mean  of  (e"rt) 

CTe-rt2  =  Variance  of  (e-rt) 

Ol2  =  Variance  of  the  cost  or  revenue  distribution. 

V 

As  previously  noted, 

Me-rt  “  E<e-rt)  =  Mt(-r), 

and 

CJ^-rt  =  M^(-2r)  -  [M^(-r)] 

Equation  4.10  can  now  be  written  asi 

V(PW )  =MC2  Mt(-2r)-(Mt(»r))  2  +  (Mt(-r))2  Cfc2 


Mt(-2r)-  (Mt(-r)| 


(4.11) 


Equation  4.11  reduces  toi 

V(Ptf)  =/Xc2[Mt(-2rHMt(-r))2]  +  ^®t(-2r)]  C^2  (4.12) 


Equation  4.9  and  4.12  are  used  to  calculate  the  mean 
present  worth  and  variance  of  present  worth  for  the 
model  elements. 

As  an  example,  let  one  assume  that  the  cost  of  a 
major  breakdown  is  normally  distributed  with  a  mean  of 
$10,000  and  a  standard  deviation  of  $100.  Furthermore,  the 
time  between  major  breakdowns  is  poisson  distributed  with 
a  mean  of  1  year.  The  nominal  interest  rate  is  10#. 
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The  expected  present  worth  is 
E(PW)  =Mc*Mt(“r) 

=  lO.OOO^e1^"* 1°“1) 

*  $9092.25  . 

The  variance  of  the  present  worth  is  calculated  from 
Equation  4.10  as 

V(PW)  =  ($10,000)2  [e(e"*20_1)-  (e(e">10-l))2] 

+  [e^e“'20_1)](100)2 
=  $760,386.05  , 

and 

0 pw  =  $872.00  . 

The  expected  present  worth  has  not  changed  from  the 
CidTis  indePendervt  modelj  however,  the  variance  has  increased 
due  to  the  cash  flow  being  a  random  variable. 

CssTss  “  Independent  Model  Element 

When  the  independent  parameters  of  cash  flow  and  time 

are  both  series  and  stochastic,  the  expected  value  of  the 

present  worth  is  derived  fromj 

PW  =  Cle-rtl  +  C2e-r(tl«2)  +  0  -r(t1+ti!+t3)  +  ____ 

°° 

E(PW)  =  E(C )  •  E  fE(e-rt)l  J 

J-l 

In  terms  of  moment  generating  functions,  the  expected 

present  worth  can  be  written  as 

E(PW)  «  |d  Mc(s)  \j  M^(-r)  (4.13) 

\  d  s  s«0/ [  . 
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Since  the  mean  cash  flow  element  would  normally  be 
given  rather  than  its  corresponding  moment  generating 
function, 


e(pw)  =/i  [jMzlL 

Ll-Mt(-r) J 


(4.14) 


The  derivation  of  the  variance  for  the  present  worth  is 
included  in  Appendix  1(4)  and  ist 


V(PW)  =  Mt(-2r)  +llc2 

Mt(-2r)- 

Mt(-r) 

2 

l-Mt(-2r) 

l-Mt(-2r)j 

l-Mt(-3 

2 

Equation  4.14  and  4.15  are  used  to  calculate  the 
expected  present  worth  and  variance  of  present  worth  for 
the  CSSTSS  model  elements. 

As  an  example  problem,  let  us  assume  that  the  damage 
to  a  bridge  during  a  flood  season  is  exponentially  dis¬ 
tributed  with  a  mean  of  $2,800,000  and  variance  of  $200,000. 
The  time  between  major  floods,  which  cause  this  damage,  is 
poisson  distributed  with  a  mean  of  20  years  and  variance  of 
10  years.  The  nominal  interest  rate  is  25$. 

By  Equation  4.14  the  expected  present  worth  is  cal¬ 
culated  as « 


E(PW)  mjl  Mt-(r) 


l-Mt(-r) 


=  2,800,000  £ 


20(e-25-l) 


l-e20(e-'25-l) 


15) 


=  $33,969*03 
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and  the  variance  is  calculated  using  Equation  4.15  as« 
V(PW)  =  $1,932,278,896. 


The  standard  deviation  is  $43,957.71.  Using  totally 
deterministic  data,  the  expected  present  worth  could  be 
calculated  fronu 


PW 


00 


-  Me  I  [e-r<20)]n=  fl 


n=l 


p-20r 


=  2,800,000 


e-20(.25) 

1_e-20(.25) 


=  2,800,000  (.00678) 

=  $18,994.23  e 


The  error  in  the  expected  present  worth  ist 


33.969.03  -  18,994.23 
33.969.03 


x  100  =  kUr.1% 


GisTcd  ”  Independent  Model  Element 

The  cash  flow  model  elements  which  are  impulsive  and 
stochastic,  with  time  being  continuous  and  deterministic, 
can  be  found  in  examples  of  monthly  receipts  or  revenues 
where  the  cash  flow  per  month  is  expressed  as  a  random 
variable.  The  expected  present  worth  is  developed  by  once 


again  looking  ati 

PW  =  Ce~rt 
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Since  t,  in  this  case,  is  known j  then 

E(PW)  =  E(C )  e’rt  «  Mce“rt  (4.16) 

One  should  note  that  this  is  the  first  model  developed 
where  common  techniques  of  using  the  mean  of  a  random 
variable  in  a  deterministic  manner  yields  a  correct 
solution. 

Since  t  is  a  constant,  the  variance  for  the  present 
worth  is  expressed  ast 

V(PW)  =  Oq2  e~2rt  . 

CssTCd  -  Independent  Model  Element 

When  the  cash  flow  is  series,  stochastic  and  the  time 
is  continuous,  deterministic,  then  the  present  worth  is 
expressed  asi 

PW  =  Cie-rtl  +  c2e-r(tl+t2)  +  C3e’r(tl+t2+t3) 

+  C4e-r^tl+t2+t3+t4)  +  .  (4.17) 

If  one  considers  that  each  time  interval  is  equal,  then 

Equation  4.17  reduces  to 
co 

Z  C i  (e-rt)j  . 


PW  = 


(4.18) 
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When  one  takes  the  expected  value  of  both  sides  of  Equation 

4.18,  while  realizing  that  E(Cj)  =  E(Ck)  =  E(C)  for  all  j 

and  k,  then  one  has» 

oo 

E(PW)  =  Y*  E(C)  (e“rt)j  , 

3*1 

where  e  is  a  constant,  the  expected  present  worth  can 
then  be  expressed  as» 

E(PW)  =  fJL  -e~r-^  .  (4.19) 

c  l-e“rt 

Letting  t  =  1,  then 

E(PU)  =  llc 

1-e  L 
or 

E(PW)  = 

i 

which  is  the  traditional  expression  for  the  present  worth 
of  an  infinite  stream  of  equal  end-of-period  cash  flows  at 
effective  interest  i. 

If  the  time  periods  are  not  equal,  then 

&  f 

E(PW)  «  fl  L  e"ri=lti  , 

0  3=1 

and  this  calculation  can  be  accomplished  by  finding  the 
mean  and  variance  of  the  distribution  of  times  and  using 
equations  for  the  CssTsg  model  elements.  Since  the  model 
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would  then  be  a  CSSTSS  model,  the  CggTC(j  model  must  assume 
constant,  equal  time  periods. 

The  variance  for  the  present  worth  of  a  cssTcd  model 
element  can  be  derived  in  a  lengthy  manner  as  was  accomplished 
for  the  CggTgg  variance,  or  one  can  simply  note  that  the 
moment  generating  function  for  a  constant  is 
M^(-r)  =  e-1^ 


Using  the  expression  for  variance  of  present  worth  in 
the  css:rss  model  (Equation  4.15),  and  substituting 

Mt(-2r)  =  e"2rt  and  M^(-r)  =  e“r”t  , 
leads  to  the  expression  for  the  variance  of  the  present 
worth  for  the  CsgTcd  model  element ; 


!_e-2rt  (a_e-2rt)(i_e-rt)2 


(4.20) 


Equation  4.20  reduces  to; 

V(PW)  =  O'  2  JSZ 

c  1-e 


2rt 


(4.21) 


An  example  of  the  CggTccj  model  element  is  dispersments 
for  operating  expenses  which  are  received  at  random  points 
during  a  month  and  paid  at  the  end  of  each  month. 


^is^cs  ~  Independent  Model  Element 


When  the  cash  flow  elements  are  impulsive  and  series  in 


nature  with  the  time  between  occurrences  being  continuous 


and  stochastic,  one  should  immediately  note  the  parallel 
with  the  model.  The  only  difference  is  time  being 

treated  as  a  continuous  element.  When  dealing  with  moment 
generating  functions,  the  distinction  between  discrete  and 
continuous  time  is  not  required  in  the  model.  The  expected 
present  worth  and  variance  of  present  worth  are  (from 
Equations  4.9  and  4.12): 

E(PW)  =  Mc Mt(-r) 

and 

V(PW)  =  jJLcZ  [Mt(-2r)  -  (Mt(-r))2]  +  Mt(-2r)  •  O'  2 

respectively.  Table  4,3  can  be  used  for  the  moment  generating 
function  for  the  time  variables  by  simply  substituting 
M^-r)  for  Mc(-r)  in  the  table.  When  the  moment  generating 
functions  for  a  specific  time  distribution  are  not  available, 
then  one  must  calculate  the  expected  present  worth  as: 

E(PW)  =  /ic  Je~rt  f(t)dt  . 

The  summation  sign  in  the  CigTis  model  has  been  replaced 
with  an  integral  sign  since  one  now  has  a  continuous 
distribution  to  deal  with.  A  similar  substitution  is 
required  to  calculate  variance.  Once  again,  it  should  be 
pointed  out  that  this  research  assumes  the  distributions 
of  interest  (both  continuous  and  discrete)  have  moment 
generating  functions. 
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CSgTcs  -  Independent  Model  Element 

When  cash  flow  is  series  and  stochastic  and  time  is 
continuous  and  stochastic,  the  model  is  the  same  as  the  one 
developed  for  CSSTSS  and  from  Equations  4.14  and  4.15« 

e(pw)  =  fi  ,M-rL  , 
c  l-M^-r) 


V(PW)  =  6  2  rM-2r)  +  f,  2 

Mt(-2r)- 

Mt(-r)i 

1  2 

c  — .  .  r'-c 

l-Mt(-2r) 

L.  | 

l-Mt(-2r) 

l-(Mf 

-r) 

|2 

Once  again,  it  is  assumed  that  the  moment  generating 
functions  for  the  continuous  time  distributions  exist. 


Summary  of  Independent  Model  Elements 

Table  4.4  is  a  summary  of  the  independent  model 
elements.  Although  each  model,  has  been  derived  separately 
with  infinite  cash  flow  streams,  one  should  note  that  all 
models  are  special  cases  of  one  of  two  general  models. 

The  CidTis,  CisTC(i,  and  CisTcs  are  all  special  cases  of  the 
Cj.gT^g  model.  Also  Cg^jTgg,  CggTcd,  and  CggTgg  are  special 
cases  of  the  CsgTss  model.  The  CssTgs  model  is  extended 
in  the  next  section  for  cases  where  cash  flow  streams  are 
finite . 
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Extension  for  Finite  Cash  Flow  Streams 

The  series,  independent  model  elements  assume  an 
infinite  number  of  cash  flows.  Although  these  formulas  may¬ 
be  fairly  accurate  when  performing  analysis  on  long-life 
items  and/or  analysis  with  high  interest  rates,  these 
formulas  would  not  be  valid  for  analysis  on  shorter-life 
items  and/or  analysis  with  low  interest  rates. 

Rather  than  developing  four  separate  models  for  finite 
cash  flow  streams  for  the  CsdTss,  CssTcd,  CgsTcs,  and 

cssTss  model  elements,  only  the  CgsTss  model  element  for 

finite  cash  flows  is  developed.  This  is  due  to  the  fact 
that  the  remaining  three  model  elements  are  special  cases 
of  this  model  element. 

The  present  worth  for  the  CgsTsS  finite  cash  flow  model 


elements  can  be  expressed  as: 

PW  =  +  C2  e-r(ti+t2) 

+  c3  e-r^l+t2+t3>  + 


n 


0 

y  l-rt 

=  C  4  e1-* 


3-1 


+  Cn  e-r(tl+tZft3+’**tn 


Since  both  the  C^’s  and  t^'s  are  each  independent,  identically 
distributed  random  variables, 


E(PW)  *  E(C ) 


£  [E(e-rt)]k 
k=l 
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Letting  X  =  E(e”rt),  the  summation  term  can  be  expressed 
ast 


n 

sn  =  Z  xk  =  x  +  X2  +  X3  +  . . .  Xn  . 


k=l 


Multiply  Sn  by  X  and  then  subtracting  from  Sn  leads  to 
Sn(l-X)  =  X-X1^1  . 

Hence , 


e;  _  X-Xn+1 
£>n  -  - 

1-X 


Substituting  X  =  E(e”rt)  leads  to 

.  E(e-rt)-|E(e-rt)|rH-1 


E(PW)  =  E(C) 
Using  S(C)  =  U 


1-E(e'rt) 


and 


E(e~rt)  =  Mt(-r) 


leads  to  the  expected  present  worth  for  the  CggTss  finite 


cash  flow  stream  model  element i 

Mt(-r)-(Mt(-r)j  n+1 


E(PW)  =  U 


(4.22) 


l-Mt(-r) 

To  develop  the  variance  of  present  worth  for  the  model • 
one  must  first  find  an  expression  for  E(PW^). 
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(PW2)  =  Cl2  e-2rt!  +  ^  e-2r(ti+t2) 

+  C32  e-2r(ti+t2+t3)  +  ...  +  Cn2  e-2r(ti+t2+t3+...tn) 

+  20^2  e-2rtl-r't2  +  20^3  e-2rtl-rt2-rt3  +  ... 

+  2C2C3  e-2rtl-2rt2-rt3  +  ... 

+  2Cn-1Gn  e-2rtl-2rt2-2rt3 - 2rtn-l-rtn  # 

Taking  the  expected  value  of  (PV/2)  and  usingi 
E(Ci*Cj)  =  E(C)  2  for  all  i  and  j , 

E(e“rti)  =  E(e“rt0)  =  E(e"rt)  =  Y, 

and 

E(e*2rti)  =  E(e“2rtj)  =  E(e_2rt)  =  X, 
leads  to 

n  n-1  n-k 

E(PW2)  =  E(C  j2)  ^  X**  +  2 [e(C)].2  Y  Y,  xkYj* 

0=1  k=l  j=l 

Using  V(C)  =  E(G2)-  [E(C)]  2, 


leads  to 

E(PW2)  = 


n  n-1  n-k 

v(c)+  |e(c)J  2J  Y  x^  +  2 (e(c ))  ZY  Y  xkY-*  * 

J  0=1  k=l  j=l 


Now  using 


n 

sn  =  Y  XJ  =  X  +  X2  +  X3  + 

j=l 


•  •  • 


+  X 


n 


mHMMUM 
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Sn(l-X)  =  X  -  X 


n+1 


s  =  X  -  x 
2>n  - 

1  -  X 


n+1 


and 


Gn  =  t  L  *3  "  X'Yl  +  ^  +  xV  +,,‘+  X  Y 

k=l  j=l 

+  x2y*  +  x2y2  +  x^y3  +...+  x^ 


l„n-l 
n-2 


+  X3Y^  +  X3*2  +  X^3  +•••+  x^  ^ 


;  xn_2Y  +  xn-2Y2 


+  Xn-*Y 

X  r  _  y2  x^y  +  x2Y2  +•  •  •+  X^Y 
-  un 

Y  ,  , 

+  X3  +  x3y  +  x3y2  +...+  X^H'3 


+  x4  +  X^Y  +  xV  +...+  X  Y 
+ 

;  xn_2  +  xn_2x  +  x"-^2 
+  x11-1  +  xn_1x 


4.n-4 


+  X 


n 


* 


Gn(£  -  1)  =  X2  +  +  X4  +...+  Xn 

-  [xV  +  xh2  +  xV  +...+  xV1-1 

n  n-1 

=  £>3  -  x  I*3- 

j=2  j=l 


n  n-1 


Y 


leads  to 


E(PW2) 


(V(c)  +  [e(C)]2)*( 


) 


Now, 

=  [Mt(-2r)]Z  -  [Mt(-2r)]n+1 
j=2  1  -  Mt(-2r) 


X  ^yj  _  M^.(-2r)  M^-(-r)  [l-  M^-(-r)] 
j=l  1  -  Mt(-r) 


and 
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Expressions  for  finite  cash  flow  streams  for  the  ramp, 
decay,  and  growth  are  developed  later  in  this  chapter. 

Dependent  Model  Elements 

Cash  flow  elements  in  this  .section  assume  dependence 
with  time.  The  common  ramp,  decay,  and  growth  functions  are 
used  as  dependence  models.  These  functions  are  widely  used 
in  economic  analysis  (ramp  functions  for  maintenance  and 
deterioration,  decay  functions  for  startup  and  learning 
costs,  and  growth  functions  for  wearin  maintenance  costs). 
Figure  4.1,  4.2,  and  4.3  illustrate  sample  ramp,  decay, 
and  growth  functions. 


Time 


Figure  4.1 
Ramp  Function 


^.idTis  “  Dependent  Model  Elements 

The  models  for  dependence  between  cash  flow  and  time 
where  cash  flow  is  impulsive,  deterministic  and  time  is 
impulsive,  stochastic  are  presented  for  ramp,  decay,  and 
growth  functions. 

C'id'^is  -  Dependent.  Ramp  Model  Element.  The  present 
worth  for  this  model  is« 

PW  =  Ct  e"rt  (4.24) 

The  expected  present  worth  can  be  stated  as« 

E(PW)  =  E(C)  •  E(t  e“rt ) .  (4.25) 

In  this  case,  the  cash  flow  (C)  is  deterministic,  and 
£(C)  =  C.  Also, 

* 

E(t  e~rt )  =  E  e“rt  =  -Mj.'  (-r). 

» 

Hence,  Equation  4.25  can  be  reduced  to« 

E(PW)  =  -CMt'  (-r) 
for  the  expected  present  worth. 

To  derive  the  expression  for  the  variance  of  the 
cidTiR  dependent,  ramp  model  element  one  must  first  derive 

p 

an  expression  for  E(PW  )i 

E(PW2)  =  E  (Cte-rt)2  »  E(C2t2e"2rt)  . 


(4.26) 
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However,  one  notes* 

Kt(-2r)  =  E(e'2rt), 

Mt'(-2r)  *  E(-2te’2rt), 


and 

Mt"(-2r)  =  E(4t2e"2rt). 

As  a  result.  Equation  4.26  can  be  expressed  as 

E(PW2)  =  g2Mt  ^~2r^  , 

4 


and  the  variance  for  the  model  can  then  be  expressed  as 


v(pw)  =  C2 


Mt"(-2r)  -  (Mt‘(-r)) 
4 


(4.2?) 


CidTis  -  Dependent,  Decay  Model  Element.  The  present 
worth  for  the  CidTis  dependent  decay  model  element  isi 

PW  =  G  e"at  e_rt 


=  C  e-(a+r)t  (4.28) 

Recognizing  that 

E(e-(a+r )t,  .  H^[-(a+r)|  , 

and  taking  the  expected  value  of  both  sides  of  Equation  4.26 
leads  to  the  expected  present  worth  for  the  model  element* 

E ( PW )  =  C  Mt  [  -(a+r) j  . 
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The  variance  for  the  model  element  can  be  derived  by 
noting  that 

PW2  =[c(e'^a+r^t)]  2 

E(PW2)  =  e(c2  e“2(a+r ) 

=  C2  Mt  |-2(a+r)) 

Hence, 

V(FW)  =  E(FW2)-  [e(PW)]  2 
=  C2!  Mtj-2(a+r)j 

GidTis  "  Dependent,  Growth  Model  Element.  The  present 

worth  for  the  G^dT^g  dependent  growth  model  element  is 
composed  of  two  previously  mentioned  model  elements  (step 
and  decay)* 

PW  =  C  e“rt  -  C  e"at  e"rt  . 

The  term  C  e"r"fc  is  the  step  function  and  C  e“at  e-rt  is  the 
decay  function.  When  these  two  terms  are  subtracted,  they 
become  the  growth  function* 

PW  =  C( l-e‘at)  e‘rt  . 


-  1  f—  ( a+r ))  i 


2 

j  i 


(4.29) 
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The  expected  present  worth  is  found  by  taking  the  expected 
value  of  both  sides  of  Equation  4.29  and  solving  in  terms 
of  moment  generating  functions* 

r  ^ 

E(PW)  =  C  Mt(-r)  -  Mt(-(a+r))j 
The  variance  is  calculated  as  previously  done  by  first 
calculating  E(PW2)i 

PW  =  G  e"rt  -  e^a+r)t 
PW2  =  C2  [  e"2rt  -  2e“(a+2r)t  +  e“2(a+r)t] 

E(PW2)  =  G2  Mt(-2r)  -  2Mt(-(a+2r)j  +  Mt(-2(a+r) 

V(P W)  =  E(PW2)  -  [e(PV/)]  2 

=  C2  [Mt(-2r)  -  2Mtj-(a+2r)]  +  Mt(-2(a+r) ) 

-  |t.lt(-r)  -  Mt(-(a+r)])2] 

Table  4.5  is  a  summary  of  the  C^dT^s  dependent  model 
elements. 

cisTis  "  Dependent  Model  Elements 

The  model  for  the  present  worth  of  the  G^ST^S  dependent 
model  elements  assumes  cash  flow  as  an  impulsive,  stochastic 
element  (expressed  as  a  function  of  time).  The  time  para¬ 
meter  for  the  model  elements  is  assumed  to  be  impulsive, 
stochastic.  The  means  and  variances  for  the  cash  flow  and 
time  parameters  are  considered  as  known  and  from  known 
distributions. 
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^isTis  “  Dependent,  Ramp  Model  Element.  The 

dependent,  ramp  model  element  has  a  present  worth  which  is 
expressed  as 

PW  =  Ct  e“rt  , 

where  C  and  t  are  both  random  variables. 

The  expected  present  worth  is  them 
E(PV/)  =  E(Ct  e~rt)  . 

By  the  mathematical  statement  in  Equation  4.24  for  the 
^idTis  dependent,  ramp  model  element, 

E(PW )  =  E(C)  •  E(t  e”rt) 

=  MC  .  (-Mt'(-r)) 

« 

=  -Mc  . 

To  derive  the  E(PW2)  term  in  the  variance  expression 
for  the  dependent,  ramp  model  element i 

E(PW2)  =  E  [(Ct  e“rt)2]  ,  (4.30) 

one  must  first  rewrite  Equation  4.30  as« 

E(PW2)  =  E(C2)  •  E(t2  e"2rt)  . 

By  definition  of  variance, 

0-o2  =  E(C2)  -  [e(C)]  2  . 

E(t2  .’2rt)  =  Mt"t-2r>  . 

4 


Also, 
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Hence 

E(PW2)  =  (CTc2+/Z  2)  M*  (-2r)  . 

4 

As  a  result, 

V(PW)  =  E(PW2)  -  [e(PW)]  2 

^6c^nc2)<l ±2-  [Me  «t'<-r>]2 

It 


=  Me2  V(-2r)  -  (Mt'(-r))2 


+  cf.2  . 

4 


CisTis  -  Dependent,  Decay  Model  Element.  The  present 

value  for  the  G-  T^_  dependent,  decay  model  element  is 

Xo  -IS 

expressed  ass 

PW  =  C  e~at  e~rt  , 

and  the  expected  present  worth  is  expressed  as» 

E(PW)  =  E(C )  *  E  (e~^a+r^t ) 

=  Mc  ’  Mt(-(a+r)). 

The  expression  for  E(PW2)  in  the  variance  for  the  model 
element  is 

E(PW2)  =  e(c2  e“2(a+rH 

=  E(C2)  •  E|e-2(a+r)t) 

=  (CTc2+  Mc2)  •  Mt(-2(a+r)j  . 

Hence,  the  variance  for  the  model  is 

V(PW)  =  (^c2+Mc2)  *  Mt(-2(a+r))  -  Mc2K("( a+r)))2 

=  Mc2[Mt  |-2(a+r)|  -  (Mt  (-(a+r)j  ]  2]  +  C^2^  |-2(a+r))  . 
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cisTis  ~  Dependent,  Growth  Model  Element.  The  present 

worth  for  the  C^sT^g  dependent,  growth  model  element 

(impulsive,  stochastic  cash  flow  and  impulsive,  stochastic 
timing)  can  be  represented  as« 

PW  =  C(l-e-at)  e“rt  .  (4.31) 

The  expected  present  worth  can  be  calculated  by  reformulating 
Equation  4.31j 

PW  =  C(e"rt  -  e-(a+r)t)  , 
and  then  taking  expected  values  to  yieldi 
E(PW)  =  E(C)LE(e“rt  -  e-(a+r)ljj 

•  E(C)  [  E(e_rt)  -  E(e"^a+r^t)] 

=  /lc  [Mt(-r)  -  Mt[-(a+r))]  .  (4.32) 

Using  the  standard  definition  of  variance,  one  must 
now  obtain  a  term  for  PW2, 

PW2  =  O2^-^  -  +  e-2(a+r)t]  _ 

The  expected  present  worth  is  then, 

E(PW2)  =  E(C2)  [E(e’2rt)  -  ZEfe-Ca+^'H)  +  E(e"2( a+r)t)] . 
This  expression  can  be  transformed  into 

E(PW2)  «  (^2+Mc2)  [wt(-2r)  -  2Tflt(-(a+2r))  +  Mt(-2(a+r) )]  . 
The  variance  follows  asi 

V(PW)  =  (Cf2+Mc2)[Mt(-2r)  -  2Mt(-(a+2r))  +  rvL,.(-2(a+r) )] 

-  Mc2[Mt(-r^  "  Mt(-(a+r))]2  . 

Table  4.6  is  a  summary  of  the  expected  present  worth 
and  variance  of  present  worth  for  the  C^gT^g  dependent 


model  elements. 
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C^sTCd  -  Dependent  Model  Elements 

The  presentation  for  the  CdsTcd  dependent  model  elements 

will  follow  the  same  format  as  presented  in  the  Ci<jTis 

dependent  model  elements.  In  the  CisTC(j  model,  the  time  is 
continuous,  deterministic.  Since  the  cash  flow  remains 
impulsive,  stochastic j  one  must  only  determine  the  differences 
between  the  two  models  where  the  time  parameter  is  deter¬ 
ministic  . 

CisTcd  ~  Ramp  Model  Element.  The  present 

worth  for  the  C^sTcd  dependent,  ramp  model  element  isi 
PW  =  Ct  e“rt  , 

v»4* 

where  t  e  is  a  constant. 

From  Equation  4.25  of  the  CidTis  dependent;  ramp  model 

element,  the  expected  present  worth  of  the  CisTcd  dependent , 

ramp  model  element  is 

E ( PW )  =  E(C)-(t  e'rt) 

-  Mc  t  e-rt  . 

Since  t  e_r't  is  a  constant,  the  variance  for  the 
^idTcd  dePendent,  ramp  model  element  is 
V(PW)  «  Cf02  (t  e~rt)2 

-  (T02  t2  e‘2rt  . 
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GisTcd  -  Dependent.  Decay  Model  Element.  The  present 
worth  for  the  C^sTC(j  dependent,  decay  model  element  is 
PW  =  C  e-at  e“rt 

=  c  e-<a+r)t  . 

The  expected  present  worth  is  then 
E(PW)  =  fX  •  E  e"(a+r)t 

v 

=  (JLc  •  . 

Since  e_^a+r^  is  a  constant, 

V(PW)  =  Cfc2  e~2(a+r )t  # 

2isTcd  ”  Dependent,  Growth  Model  Element.  The  present 
^orth  for  the  CigTcd  dependent,  growth  model  element  is 

PW  =  C(l-e“at)(e_rt) 

=  G(e-rt  -  e-(a+rH)  , 
where  r,  t,  and  a  are  constants. 

Since  G  is  the  only  random  variable,  the  expected 
present  worth  follows  as« 

E(PW)  =  E(0)  '  E(e'rt  -  e-(a+rH) 

=  Mc  •  [e-rt  -  e-(a+r)t]  . 

The  variance  for  the  model  is 

V(PW)  =  6QZ  [e"rt  -  e-^a+r^t]2 

=  Cfc2  [e-2rt  -  2e'(a+2r)t+  e-2(a+r)t  ]  # 
Table  4.7  is  a  summary  of  the  expected  present  worth 
and  variance  of  present  worth  for  the  C^gTcd  dependent 


model  element 
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Table  4.? 


^isTcd  DePenden,t  Model 

Elements 

Model 

E(Ptf) 

V(PW) 

Ramp 

Mo  » 

O'  2  t2  e”2rt 

C 

Decay 

Mc  e-(a+r)t 

Cfc2  e-2<a+r)t 

Growth 

Mc'[e- rt  _  e-(a+r)t] 

0'c2[e“2rt  -  2e‘(a+2r)t 

+  e-2(a+r)t] 
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CisTcs  -  Dependent  Model  Elements 

As  previously  noted  in  this  research,  treating  time  as 
continuous  or  discrete  does  not  alter  the  mathematical 
development.  This  was  due  to  the  use  of  moment  generating 
functions  which  are  assumed  to  be  known  for  both  discrete 
and  continuous  distributions.  The  expected  present  worth 
and  variance  of  present  worth  for  each  of  the  C^sTcg 

dependent  model  elements  will  be  the  same  as  for  the  Cj_sT^s 
dependent  model  elements. 


cisTcs  "  Dependent,  Ramp  Model  Element.  The  expected 


present  worth  is 

E(PW)  =  -jlc  Mt'(-r)  , 

and  the  variance  of  present  worth  is 

Mt"(-2r) 


v(? w)  -  M<-2 


-  [Mt’(-r>] 


+  6. 


2 


GigTCs  “  Dependent,  Decay  Model  Element.  The  expected 

present  worth  is 

E(PW)  =  •  Mt(-(a+r))  , 

and  the  variance  of  present  worth  is 

V(PW)  =  Mc2  [  Mt(-2(a+r))-[Mt(-(a+r))]2] 

+  0 c2  Mt(-2(a+r))  . 
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Cis'Ic  s  “  Dependent,  Growth  Model  Element.  The  expected 
present  worth  is 

E(PW>  =  Mc  K<-r)  -  M,  (-(a+r)  )  ]  , 

and  the  variance  of  present  worth  is 

V(PW)  =  (  0'c2+Mc2)[r>It(-2r)  -  2Mt(-(a+r))+  Mt(-2(a+r))] 

-  Mc2  [ Mt(-r)  -  Mt(-(a+r))]2  . 

For  consistency  in  the  presentation  of  models,  Table  4.8 
contains  the  CdsTcs  dependent  model  elements. 

Summary  of  Impulsive,  Dependent  Cash  Flow  Elements 

The  C^T^g.  CdgTcd,  and  C^sTcs  models  presented  in  this 

research  under  the  dependent  model  elements  are  all  special 
cases  of  the  CdgTis  dependent  model  elements.  This  con¬ 
clusion  parallels  the  findings  for  the  independent  model 
elements. 

The  expected  present  worth  will  be  derived  for  all  of 
the  dependent,  series  models  in  the  next  section.  The 
variance  will  be  derived  for  the  csdTss  dependent ,  ramp 

i 

function. 


CsdTss  -  Dependent  Model  Elements 

The  models  for  dependence  between  cash  flow  and  time 
when  the  cash  flow  is  series,  deterministic  and  the  time 
is  series,  stochastic  are  presented  in  this  section  for 
ramp,  decay,  and  growth  functions. 

^sdTss  -  Dependent,  Ramp  Model  Element.  The  present 
worth  for  this  model  is 

PW  =  C-ti  e~r\  +  C(t1+t2)  e"r^tl+t2^ 

+  ...+  C(t1+t2+. .  .+tn)  e-r(tl+t2+#  *  ,+tn) 

•  y 

where 

t ^  is  a  random  variable, 

G  is  a  constant. 

The  present  worth  is  reformulated  asi 

PW  .  Ct,  e-rt.[i  +  e'rt2  +  e-r(t2+t3)  +  e-r(t2+t3+t4)  +  _ 

+  Ct2  e-rt2  e-rt' [1+e-rt3+e-r(t3+t4)+e-r(t2+t3+t4)+.-,' 

+  Ct3e-rt3e-rt2e-rtl  [l+e'rtlt+e-r<t'*+t5)+e-r(t,‘+t5+t6)+...] 

•  •  •  i 

The  expected  present  worth  can  now  be  written  as 

oo  oo 

E(PW)  =  C  E(te“rt)  Y.  [E(e"rt)]n*Z  [E(e~rt)]n 

n=0  n=0 


Ill 


This  is  proven  by  noting 

E(e‘rti)  =  E(e"rtJ)  for  all  i,j 

and 

E(tie'rti)  =  E(tje"rtj)  for  all  i,j. 

Analytically,  Equation  4.33  is  not  too  difficult  to 
use?  however,  the  expression  for  the  variance  of  present 
worth  for  the  csdTss  dependent,  ramp  model  element  is  a 
great  deal  more  detailed.  As  before,  to  calculate  the 
variance,  one  must  first  generate  PW2  and  then  calculate 
E(PW2).  This  calculation  is  contained  in  Appendix  I  (5) 
and  the  resulting  variance  for  the  model  is« 

r  M  "(-2r)  M  ' ( - 2r ) 2  ir  2K+(-r)  1 

l.  4[l-Mt(-2r)]  2  [l^Mt(-2r>]3  JL  l-Ht(..r)J 

+  Kt'(-2r)Mt(-r)  _  [Mt’(-r)]2 

"~[l-Mt(-2r)]  [l-Mt(-r)]  3  ”[l-Mt(-r)]4  J.  (4-34) 

Although  this  variance  expression  can  be  further 
simplified,  it  is  complete  enough  at  tnis  point  to  draw  a 
conclusion  as  to  the  usefulness  of  such  an  expression.  The 
original  intent  of  the  analytical  solution  approach  was  to 
develop  models  which  were  mathematically  correct  and 
realistic  to  use.  Without  the  use  of  simulation,  it  would 
appear  that  calculating  variance  for  dependent,  series 


V(PW)  =  C2 
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models  is  not  of  significant  value.  The  remainder  of  the 
GsdTss  dePendent  models  and  the  remainder  of  the  series, 

dependent  models  will  be  limited  to  the  development  of 
expected  present  worths. 

C  .T  -  Dependent,  Decay  Model  Element.  The  present 

SQL  SS 

worth  for  the  CsdTgs  dependent,  decay  model  element  can  be 
represented  ass 

PW  =  Ce-<»+r>tl  +  0e-<a+r)(tl+t2)  +  Ce-(a+r)(t1+t2+t3>+_ _ _ 
Noting  that 

E(e-(a+r)ti)  =  E(e-(a+r)tj)  for  ^  k  ^ 

and 

e-(a+r)ti  independent  of  for  all  i 

and  j, 

one  can  express 

Ere-(a+r)(ti+t2+...tn)] 

as 

E[e-(a+r)t  ] n. 

Hence i 

oo 

E(PW)  =  I  o[E(e-(a+r)t)]n 
n»l 

E  (e-<a+rH) 

=  °  l-E(s-(a+r)t) 

M+(-(a+r)) 

SS  C  *  * 

l-Mt(-(a+r)) 


I 
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C'sdT’ss  -  Dependent,  Growth  Model  Element.  The  present 


worth  for  the  CsdTgs  dependent,  growth  model  element  can  be 
represented  as» 

PW  =  c[e-rt  -  e*(a+r)t]  +  c[e“r{tl+t2)  -  e-U+r)(t1+t2)] 

+  C[e'r('fcl+t2+t3)  -  e-(a+r)(tl+t2+t3}]  +  ...  . 

As  pointed  out  earlier  in  the  research,  the  growth  function 

is  the  sum  of  a  step  function  and  a  decay  function.  As  a 

result,  the  expected  present  worth  can  be  expressed  asi 
co  oo 

E(PW;  =  C  E  [E(e~rt)]n  -  C  I  [E(e"(a+r) )] n 

n=l  n=l 

Mt(-r)  Kt(-(a+r)) 

[l-M^-r)  l-Mt(-(a+r))  J  . 

The  expected  present  worth  for  the  CgdTgg,  dependent 
model  elements  are  presented  in  Table  4.9. 

CgsTgg  -  Dependent  Model  Elements 

The  models  for  dependence  between  cash  flow  and  time 
when  the  cash  flow  is  series,  stochastic  and  the  timing  of 
that  cash  flow  is  series,  stochastic  are  presented  in  this 
section.  Only  the  expected  present  worth  for  the  ramp, 
decay,  and  growth  models  are  derived. 


Table  4.9 

csdTss  Dependent  Model  Elements 


Model 


Ramp 


Decay 


E(PW) 

-  C  M^' (-r) 

(l~Mt(-r))2 

C  Mt(-(a+r) 
l-Mt(-(a+r) ) 

Mt(-r)  Mt(-(a+r)) 

.  l-Mt(-r)  l-Mt(-(a+r)) 


Growth 


C 
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CSgTss  -  Dependent,  Ramp  Model  Element.  The  present 


worth  for  the  CssTgg  dependent,  ramp  model  element  can  be 
represented  as 

PW  =  C1(t1e“rtl)  +  C2(t1+t2)e’r^tl+t2) 

+  C3(t1+t2+t3)e‘r(tlH-b2+t3)  +  ...  .  (4.35) 

The  expected  present  worth  is  derived  in  Appendix  I  (6) 


and  is  shown  to  bei 


E(PW)  =  M, 


-Mt  (-r) 

[l-Kt(-r)]2  J  . 


C-.T..  -  Dependent,  Decay  Model  Element.  The  present 

S  S  ob 

worth  for  the  CSSTSS  dependent,  decay  model  element  is 
PW  =  Oj  G-(a+r}ti  +  e-(a+r)(t1+t2)  +  ... 

+  Cn  e-(a+r)(tl+t2+-"+tn)  +  ...  . 

Since  the  t^'s  and  G^'s  are  independent  and  identically 
distributed  random  variables, 

E(PW)  •  fJL c  [  [E(e-(a+r)t)]n 

n=l 

^  Mt(-(a+r)) 

C  l-Mt(-(a+r))  . 


CSSTSS  -  Dependent,  Growth  Model  Element.  The  present 


worth  for  the  CsgTss  dependent,  growth  model  element  is» 

PW  =  Cjd-e^lKe-1^)  +  C2(l-e-a(tl+t2))e‘r(tl+t2)  +  ... 

+  Cn(l-e"a(tl+t2+,,,+t^)e-r(tl+t2+*‘*+tn)  +  ... 

=  Cl(e-rtl-e-(a+r)tl)  +  c2(e'r(tl+t2)-e"(a+r)(tl+t2)) 

+  ...  +  Cn(e"r^tl+t2+*“+tn)_e-(a+r^fcl+t2+***+tn^)  +  ... 

Taking  expected  values,  while  noting  that 
E(Ci)=  E(C  j )  =  yUc 


and 


E(e”rti)  =  E(e"rtj)  =  M^(-r)  for  all  i  and  j  , 
leads  to 


E(PW)  =  flc 
-  Me 


OO  CO 

Z  (Mt(-r))n  -  Z  [(Mt(-(a+r))]n 


n=l  v  n=l 

Mt(-r)  _  Mt(-(a+r)) 


Ll-Mt(-r)  l-Mt(-(a+r)) J  . 

The  expected  present  worth  for  the  CssTgs,  dependent  model 
elements  are  presented  in  Table  4.10. 


CSgTCd  -  Dependent  Model  Elements 

The  models  for  dependence  between  cash  flow  and  time 
when  the  cash  flow  is  series,  stochastic  and  the  time  is 
continuous i  deterministic  are  presented  in  this  section 
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for  ramp,  decay,  and  growth  functions.  Working  from  the 
previous  C__T__,  dependent  model  elements,  one  need  not 

oo  So 

derive  these  models  directly. 


CssTCd  -  Dependent,  Ramp  Model  Element.  The  CssTC(j 


dependent,  ramp  model  element  is  indirectly  derived  from 
the  C  T  dependent,  ramp  model  element.  Recalling  that 

OO  SS 

the  expected  present  worth  for  the  CggTgs  dependent,  ramp 
model  element  is 

"-M+* (-r) 

E(PW)  =  IX  -1 - , 

c  L  [l-'H(-r)]  2j 

one  need  only  recognize  that  constants  do  have  moment 
generating  functions  (e.g.  for  time  as  a  constant) i 
Mt(-r)  =  e“rt  , 

-Mt'(-r)  =  t  e"rt  . 

Hence,  the  expected  present  worth  for  the  CggTccj  dependent, 
ramp  model  element  is 


E(PW)  = 


IX cu  B'rt) 

(l-e-rt)2 


CggTcd  -  Dependent,  Decay  Model  Element.  The  CggTC(j 

dependent,  decay  model  element  is  also  derived  from  the 
CggTss  dependent,  decay  model  element. 
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^ssTcd  “  Dependent,  Growth  Model  Element.  The  CssTC(j 
dependent,  growth  model  element  is  derived  from  the  CSSTSS 
dependent,  growth  model  element  as  follows: 

E(PW)  = 

=  Me 

The  expected  present  worth  for  the  CssTc^  dependent  model 
elements  are  presented  in  Table  4.11. 

CssTcs  -  Dependent  Model  Elements 

When  cash  flo w  is  series  and  stochastic  and  time  is 
continuous  and  stochastic,  the  models  are  the  same  as  those 
presented  under  the  CSSTSS  dependent  model  elements  (see 
Table  4.10). 


Mt(-r)  Mt(-(a+r)) 
l-Mt(-r)  l-Kt(-(a+r) ) 
„-rt  -(a+r)t  ] 


.  -  ( a+r )  t 


120 


Model 


Ramp 


Decay 


Growth 


Table  4.11 

CssTC{j  Dependent  Model  Elements 


E(PW) 


t  e 


-rt 


(1-t  e-rt) 


2 


-(a+r)t 

U  £ - 

1-e-(a+r)t 

e-rt  e-(a+r)t 
0  l-e'rt  '  i-e-<a+rH. 
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Summary  of  Series  Dependent  Model  Elements 

The  series,  dependent  model  elements  axe  all  special 
cases  of  the  CSSTSS  dependent  model  elements  presented  in 

Table  4.10.  The  CSSTSS  dependent  model  elements  are  extended 
in  the  next  section  for  cases  where  cash  flow  streams  are 
finite. 

Extension  of  CgsTss  Dependent  Model  Elements 
for  Finite  Gash  Flow  Streams 

The  CssTsg  dependent  model  element  for  finite  cash 
flow  streams  are  presented  in  Table  4.12.  The  derivation 
of  the  models  are  found  in  Appendix  I  (7),  (8),  and  (9)* 

Expansion  of  General  Model 

The  general  model  for  a  present  value  analysis  was 
presented  in  Chapter  3  as« 

E(PV )  =  E(P)  +  E(E)  -  E(R)  -  E(S)  ,  (4*36) 

where 

E  =  Expected  value  operator 
P  =  Purchase  price 
E  =  Expense 
R  =  Revenue 
S  =  Salvage. 


Table  4.12 


CSSTSS  Dependent,  Finite  n  Cash  Flow 


Model  Elements 


Expected  Present  Worth 


Growth 
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Let  us  now  assume  that  there  are  i  independent  alter¬ 
natives  and  each  of  the  alternatives  can  have  up  to  k  model 
elements  associated  with  P,  E,  R,  and  S.  The  expected 

xy. 

present  worth  of  the  ix  alternative  is  then: 

k  k 

ECPVi)  =  Z  E(Pk)  +  Z  E(Ek) 

3*1  0=1 

k  k 

-  Z  E(R.  )  -  Z  E(Sk) 

3=1  k  3=1  k  • 

The  analyst  must  now  identify  the  type  of  cash  flow 
model  element  and  compute  the  expected  present  worth  for 
that  model  element  using  Tables  4.9  through  4.12.  Once 
this  has  been  accomplished,  the  expected  present  worth 
would  be  added  to  the  appropriate  category  of  purchase  price, 
expense,  revenue,  or  salvage  value.  Equation  4.36  can  then 
be  used  to  calculate  the  expected  present  value  for  that 
alternative. 

The  variance  can  be  calculated  in  a  similiar  manner 
and  is  equal  to 

k  k  k  k 

V(Wt)  =  i  V(Pi>  +  LV(Ek)  t  Zv(Rk)  +  £v( sk)  • 

3=1  3=1  3=1  3=1 

Many  additional  model  element  types,  not  discussed  in 
this  research,  can  also  be  included  in  the  general  model 
for  the  calculation  of  expected  present  value.  The  general 
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model  for  the  variance  of  present  value  would  also  be  valid 
as  long  as  the  cash  flows  were  considered  to  be  independent. 
Dependent  cash  flow  elements  could  be  handled  with  an 
expanded  model  as  long  as  time  is  not  a  random  variable. 

If  time  is  a  random  variable,  it  is  not  worthwhile  to  develop 
a  covariance  matrix  for  the  cash  flows  because  of  the 
difficulty  involved. 

The  actual  selection  among  alternatives  would  normally 
be  based  upon  minimum  expected  present  value.  There  are 
two  possible  cases  where  variance  of  present  value  would 
enter  the  selection  procedure.  First,  the  expected  present 
values  of  two  alternatives  could  be  either  tied  or  too  close 
for  the  analyst  to  be  positive  which  was  the  better  choice. 

In  this  case,  the  alternative  with  the  smaller  variance 
would  be  selected.  The  second  case  is  where  the  variance  of 
the  cheapest  alternative  is  greater  than  the  variance  of  an 
alternative  close  to  the  expected  present  value  of  the 
cheapest  alternative.  Although  more  research  is  required 
to  justify  this  viewpoint,  there  is  some  intuitive  appeal 
to  minimize  uncertainty  even  if  a  slightly  higher  cost  is 
involved. 

The  next  chapter  contains  an  indepth  discussion  of  the 
contributions  of  this  research. 


Chapter  5 


DISCUSSION  OF  MODELS 

This  research  has  extended  the  present  knowledge  of 
capital  investment  problems  in  several  ways.  A  detailed 
taxonomic  structure  has  been  developed  which  is  utilized  in 
this  research  as  an  outline  for  model  development.  The 
taxonomic  structure  also  serves  as  a  set  of  designations 
for  new  models  which  should  aid  in  organization,  storage, 
retrieval,  and  research  of  these  models.  Many  specific 
models  have  also  been  developed.  These  models  should  be 
useful  to  the  practitioner  and  of  motivational  value  to 
structuring  realistic  classroom  presentations.  Each  of 
these  areas  are  addressed  in  this  chapter. 

Taxonomic  Structure 

There  is  a  definite  requirement  for  a  taxonomic  structure 
for  capital  investment  problems  and  economic  analysis  prob¬ 
lems.  Fleischer  and  Ward  (1977:24)  developed  such  a  taxonomic 
structure  and  their  work  is  presented  in  Table  3*7*  The 
intent  of  their  work  was  three-fold.  First,  the  taxonomic 
structure  serves  as  a  "road  map"  for  locating  and  selecting 
appropriate  models.  Second,  the  taxonomic  structure  serves 
as  a  method  of  organization,  storage,  and  retrieval  of 
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models.  Third,  the  taxonomy  permits  enumeration  of  simpler 
models.  As  part  of  their  research  effort,  a  literature 
search  was  accomplished  to  designate  which  areas  of  economic 
analysis  had  been  researched  and  to  determine  where  addition¬ 
al  research  was  required.  Table  3.7  is  reproduced  in  Table 
5*1  with  additional  entries  for  the  contributions  of  this 
research.  The  entries  5a,  5t>»  5c,  and  5d  in  the  table 
highlight  where  the  research  was  accomplished. 

The  taxonomic  structure  developed  by  Fleischer  and 
Ward  (1977)  was  not  detailed  enough  to  permit  the  type  of 
enumeration  of  models  described  in  this  research.  Additional 
detail  was  added  to  facilitate  the  research.  A  review  of 
Table  3-7  and  Figure  3*H  highlights  the  additions  to  the 
taxonomic  structure.  The  cash  flow  has  been  subdivided  into 
expenses  and  revenues.  An  additional  category  has  been 
made  for  independent  and  dependent  models.  The  impulsive 
cash  flow  has  been  subdivided  into  those  types  of  cash  flows 
which  only  occur  once  (designated  as  impulsive)  and  those 
cash  flows  which  are  repeated  with  some  type  of  frequency 
(designated  as  series).  Additional  models  have  been  developed 
for  series  cash  flows,  both  independent  and  dependent,  where 
the  cash  flow  streams  are  finite. 

Independent  Models 

Independent  models  are  developed  in  Chapter  4  with  a 
summary  in  Table  4.4.  The  C^dT£s,  CdsTcd,  and  C^sTis  models 


I 


Table  5.1 

Taxonomic  Structure  and  Current 
State -of -the- Art 


Cash 

Flow 

Amount 

Timing 

How 

Cert 

When 

Cert 

Det 

Det 

Disc 

Sto 

<D 

> 

Det 

•H 

W 

Cont 

Sto 

3 

p. 

Det 

B 

M 

0 

-p 

Disc 

Sto 

C/2 

Det 

Cont 

Sto 

Det 

+> 

Disc 

Sto 

ra 

3 

0 

0 

Q 

Det 

3 

C 

Cont 

Sto 

-P 

C 

Det 

O 

O 

0 

-p 

Disc 

Sto 

co 

Det 

Cont 

Sto 

Interest  Rate 


Stepwise 
Constant  Constant 


Det  Sto  Det  Sto  Det  Sto 


1-Reisman  and  Rao  (1973) 

Cert -Certainty 

2-Ward  (1975) 

Det-Deterministic 

3-  Hillier  (1969) 

Sto-Stochastic 

4-Hillier  (1963) 

Cla-Classical 

5-Estes  (current  research) 

C  ont-C  ont inuous 

Disc -Discrete 

Source: 

G.  A.  Fleischer  and  T.  L.  Ward.  "Classification  of 
Compound  Interest  Models  in  Economic  Analysis."  Engineering 
Economist,  (Fall,  1977). 24. 


128 


are  shown  to  be  special  cases  of  the  Cj_sTj_s  model  which 

assumes  that  cash  flow  and  timing  of  cash  flow  are  impulsive 

and  stochastic  in  nature.  The  C„,T  .  C^JT  , ,  and  G  T 

sd  ss  ss  cd  ss  cs 

models  are  special  cases  of  the  CSSTSS  model  which  assumes 
that  the  cash  flows,  and  associated  timing  of  the  cash  flows, 
are  series  in  nature  with  stochastic  parameters.  Expected 
present  worth  and  variance  of  present  worth  are  derived  for 
each  model.  An  extension  for  the  CSSTSS  model  is  presented 
next  for  models  which  have  a  finite  cash  flow  stream.  The 
expected  present  worth  and  variance  of  present  worth  for 
the  finite  cash  flow  CSSTSS  model  are  presented  in  Equations 
4.22  and  4.23  respectively. 

These  models  are  of  interest  in  that  they  handle  data 
for  capital  investment  where  such  variables  as  initial  cost, 
cost  and  timing  of  breakdowns,  and  salvage  values  must  be 
considered  as  random  variables. 

Dependent  Models 

Dependent  models  are  also  developed  in  Chapter  4.  The 
assumption  made  in  these  models  is  dependence  between  cash 
flow  and  time.  The  ramp,  decay,  and  growth  functions  are 
developed  for  each  of  the  eight  basic  model  elements 
previously  derived  for  the  independent  models.  The  C^gT^s 


models  are  presented  in  Table  4.6.  The  remaining  impulsive 
models  are  all  special  cases  of  this  model. 


The  series,  dependent  models  (CsdTss,  CssTcd,  CssTcs) 
are  special  cases  of  the  CSSTSS  dependent  models  presented 
in  Table  4.10.  Variance  of  present  worth  is  developed  for 
the  CsdTss  dependent,  ramp  model.  As  can  be  seen  in 
Equation  4.34,  the  complexity  of  the  variance  expression 
makes  it  of  little  value  as  an  analytical  tool.  In  addition, 
attempts  to  derive  a  general  form  for  the  CssTss  variance 
were  not  successful.  Expressions  for  the  expected  present 
worth  for  the  CSSTSS  finite  cash  flow  dependent  model 
elements  are  presented  in  Table  4.12.  The  development  of 
dependent  series  models  is  considered  to  be  of  relevance 
due  to  the  use  of  ramp  functions  with  maintenance  and 
deterioration  costs,  decay  functions  with  start  up  and 
learning  costs,  and  growth  functions  with  wearin  maintenance 
costs. 


Summary 

A  total  of  nine  independent  models  are  developed  with 
expected  present  worth  and  variance  of  present  worth.  In 
addition,  twelve  dependent  models  are  presented  with  expected 
present  worth  and  variance  of  present  worth  for  impulsive 
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dependent  models.  For  the  series,  dependent  models,  only 
expected  present  worth  is  presented  for  eleven  models  and 
the  finite  cash  flow  stream  model.  A  total  of  36  models 
have  been  presented  in  Chapter  4  with  the  majority  having 
both  expected  present  worth  and  the  variance  of  present  worth 
developed. 

In  order  to  clarify  the  specific  areas  of  comparison 
between  the  taxonomic  structure  accomplished  by  Fleischer 
and  Ward  (1977:24)  with  this  research,  Table  5.1  should  be 
compared  to  Tables  5*2  through  5.5. 

The  next  chapter  contains  a  development  and  justification 
for  a  network  logic  to  solving  problems  which  utilize  the 
analytical  tools  developed  in  this  research. 


Chapter  6 


NETWORK  REPRESENTATION  OF  A  CAPITAL 
INVESTMENT  ALTERNATIVE 

All  textbooks  which  discuss  alternative  selection 

models  for  engineering  economics  develop  rate-of-return 

formulas.  These  formulas  are  developed  by  the  use  of  both 

mathematical  equations  and  visual  depictions  of  the  process. 

As  stated  by  A.  Alan  B.  Pritsker  ( 1979s 1), 

The  modeling  of  a  system  is  made  easier  if: 

1)  physical  laws  are  available  that  pertain  to  the 
system;  2)  a  pictorial  or  graphical  representation 
can  be  made  of  the  system;  and  3)  the  variability 
of  system  inputs,  elements  and  outputs  is  manageable. 

This  research  has  attempted  to  improve  upon  models  for 

capital  investment  problems  by  providing  useful  tools  to 

accomplish  the  modeling  of  an  economic  process.  The  research 

is  expanded  in  this  chapter  to  offer  suggestions  as  to  how 

a  pictorial  representation  can  be  made  of  the  process. 

Present  Worth  Development 

Given  a  future  sum  E,  a  continuous  interest  rate  r, 

-rt 

and  time  of  occurrence  t,  the  present  worth  PW  is  E  e  . 
Figure  6.1  represents  a  network  configuration  of  the 
problem.  The  nodes  of  the  network  can  be  used  to  generate 
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and  collect  dollar  values.  The  branches  are  used  to 
designate  time  between  nodes.  The  symbol  E  was  chosen  to 
designate  an  expense. 


Figure  6.1 

Present  Worth  of  a  Single  Cash  Flow 

Let  us  now  extend  this  pictorial  network  by  allowing 
an  initial  expense  P  occurring  at  t=0  and  a  series  of 
n  expenses  (Ej,  E2,  Ey...En)c  Figure  6.2  illustrates  this 
cash  flow. 


Figure  6 . 2 

Present  Worth  of  a  Multiple  Cash  Flow 
The  present  worth  is 
n 


I  Ei  e"rti  • 


PW  =  P  + 
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If  one  were  to  assume  that  all  of  the  expenses  and  the 
time  between  expenses  are  random  variables,  Figure  6.2  is 
still  valid.  Only  the  method  used  to  calculate  the  present 
worth  has  changed.  Obviously  for  simple  cash  flows,  such  as 
those  discussed  to  this  point,  a  picture  is  probably  not 
required  to  aid  the  analyst  in  solving  the  problem  at  hand. 

Now  let  us  expand  upon  the  nature  of  the  expenses  involved. 

The  following  example  problem  is  presented  as  a  more 
realistic  real-world  problem:  A  manufacturing  firm  has 
decided  to  expand  their  plant  to  produce  a  new  product. 

A  lathe  is  required  in  one  area  of  the  production  line. 

Several  lathes  which  have  passed  an  initial  screening  to 
insure  that  they  can  meet  or  exceed  minimum  production 
standards  have  been  suggested  and  are  now  ready  for  economic 
evaluation.  The  data  for  the  first  lathe  to  be  considered 
has  been  accumulated  from  the  producer  of  the  lathe  and 
in-house  estimates.  The  initial  cost  of  the  lathe  is  $50,000. 
The  installation  cost  is  dependent  upon  suspected  electrical 
problems  and  is  estimated  to  be  from  a  triangular  distri¬ 
bution  with  a  mean  of  $10,000  and  the  endpoints  of  $2,000 
and  $15,000.  The  salvage  value  for  the  lathe  is  normally 
distributed  with  a  mean  of  $5»000  and  a  variance  of  $750. 

This  particular  lathe  has  been  extensively  used  in  industrial 
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applications  and  Table  6.1  is  a  summary  of  the  repair  data 
provided  by  the  manufacturer.  Time  to  failure  is  normally 


distributed  with  a  mean  of  60  days  and  a  variance  of  10 
days.  The  cost  of  down-time  is  a  constant  $500  per  day. 
Monthly  operational  costs  for  preventive  maintenance  and 
electrical  consumption  is  $600.  The  compounding  period  is 
daily  and  the  rate  of  return  is  20 %.  Service  life  is  5 
years . 

Figure  6.3  represents  a  possible  network  type  con¬ 
figuration  for  the  problem.  The  nodes,  designated  as 

© 

(where  k  stands  for  the  node  number)  is  where  any  expenses 
generated  at  that  time  are  introduced.  Transactions  passing 
thru  the  network  would  collect  present  worth  for  all 
expenses  up  to  that  point  in  time.  The  symbol 

© 

is  used  as  a  decision  node  to  designate  that  some  parameters 
must  be  tested  prior  to  determination  of  branching.  The 
use  of  a  triangle  on  the  output  portion  of  a  node  is  used 
to  designate  probabilistic  branching.  The  absence  of  a 
time  specification  on  a  branch  assumes  that  zero  time  is 
associated  with  that  branch.  The  actual  present  worth  for 


I 
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Table  6.1 

Repair  Data  for  Lathe  1 


Type  of 
Failure 

Probability 

Cost  ($) 

Probability 
of  Repair 

Down-time 
( days ) 

Major 

.1 

N(1K,  100) 

•  90 

4 

Normal 

•  5 

u(ioo,  500) 

.75 

2 

Minor 

.4 

E( 100) 

.50 

1 

Legend: 

N  is  normal  distribution 
U  is  uniform  distribution 
E  is  exponential  distribution 


one  realization  of  the  network  could  be  carried  along  with 
cash  transaction  by  an  attribute. 


The  approach  suggested  to  solve  this  problem  is  by 
using  simulation?  however,  many  analytical  tools  (including 
those  introduced  in  this  research)  could  be  used  as  actual 
nodes  in  the  network. 

The  example  problem  could  be  extended  further  to 
approach  the  real-world  more  closely  by  adding  taxes, 
revenues,  depreciation,  etc.  As  the  complexity  of  the 
problem  increases,  the  need  for  a  pictorial  representation 
is  justified. 

Preliminary  research  indicates  that  QGERTS,  developed 
by  A.  Alan  B.  Pritsker  is  adaptable  to  performing  the  analysis. 
It  would  be  a  major  undertaking,  requiring  extensive  re¬ 
writing  and  additions  to  the  current  program?  however,  it 
is  suggested  in  Chapter  8  of  this  work  that  this  task  be 
undertaken.  Other  possible  simulation  programs  which 
appear  to  have  potential  are  SLAM  and  GPSS. 

Summary 

The  intent  of  this  chapter  is  not  to  design  a  simulation 
language  to  be  used  in  solving  capital  investment  problems. 

The  chapter  does  layout  a  network  procedure  to  use  in 
solving  problems.  The  next  chapter  introduces  some  ap¬ 
plications  using  the  results  of  this  research  and  the 
netv/ork  logic  will  be  used  as  an  illustrative  tool. 


Chapter  7 


APPLICATIONS  FOR  RESEARCH  MODELS 

There  are  three  basic  areas  where  this  research  is 
useful.  The  first  area  is  in  adding  more  realistic 
analytical  tools  to  those  already  used  in  traditional  or 
classical  engineering  economic  problems,  the  second  area 
is  in  the  validation  and  verification  of  simulation  models, 
and  the  third  area  is  in  the  direct  use  of  this  research 
in  the  development  of  a  simulation  language  for  engineering 
economics.  Each  of  these  areas  are  addressed  in  this 
chapter. 


Extensions  to  Traditional  Engineering 
Economic  Problems 

This  section  of  applications  is  concerned  with  adding 
this  research  to  the  set  of  analytical  tools  used  to  solve 
traditional,  engineering  economic  problems.  Areas  covered 
are  capitalized  cost  comparisons,  replacement  economy,  and 
lease  or  buy  decisions.  The  research  is  not  limited  in 
application  to  these  areas;  however,  these  areas  were 
selected  to  show  the  value  of  the  research  in  a  traditional 
type  problem. 


Capitalized-Cost  Comparisons 

Capitalized-cost  comparisons  are  present  worth  analysis 
where  the  comparison  period  is  assumed  to  be  infinitely 
long.  Taylor  ( 1964 1  95)  points  out  that  capitalized  cost 
predominated  in  the  early  days  of  railroad  construction 
and  expansion  due  to  the  viewpoint  that  roadbeds,  tracks, 
and  bridges  were  considered  to  have  the  "mathematical 
equivalent  of  perpetual  lives".  This  is  due  to  the  very 
small  error  introduced  in  using  discounting  factors  with 
infinite  lives  versus  long  life.  This  error  is  decreased 
even  more  with  higher  interest  rates.  Capitalized  costs 
comparisons  are  valid  for  many  analysis  problems  such  as 
tunnels,  dams,  aqueducts,  bridges,  and  interstate  highways. 

A  few  typical  analysis  problems  are  i3.1ustrated  in  the 
following  examples. 

Example  7 . 1  The  following  example  is  modified 
from  Taylor  (1964*  95,  ex.  7.8). 

A  dam  costing  $100,000  to  construct  will  cost  $1 5,500 
a  year  to  operate  and  maintain.  Another  design  costing 
$150,000  to  build  will  cost  $10,000  a  year  to  operate  and 
maintain.  Both  installations  are  felt  to  be  permanent. 

The  annual  dispersments  are  assumed  to  be  lump  sum,  year 
end  payments  and  the  interest  is  continuous  at  25#* 


1 45 


Taylor's  solution  would  be  to  calculate  the  present 
worth  of  the  first  design  ast 


$100,000  +  $15,500  * _ _  =  $154,572.58  , 

e  -°-l 


and  the  second  design  asi 

$120,000  +  $10,000  *  — 1— _ 

e  -25-i 


$155,208.12 


The  first  design  is,  in  terms  of  present  dollars,  $635*54 
better.  Now,  let  one  assume  that  the  annual  dispersments 
are  dependent  upon  funding  by  the  government  and  that  this 
funding  is  normally  distributed  with  a  mean  of  one  year 
and  a  variance  of  6  months.  The  first  cost  of  both  designs 
is  assumed  to  be,  as  before,  a  constant. 

Figure  7 .1  illustrates  a  network  for  the  decision 
process.  The  Cs(jTsg  model  element  is  used  to  calculate 
the  present  worth  of  the  first  design  asj 


$100,000  +  $15,500 


Mt(-.25) 
,1-Mt(-.25)J 


M^(-.25)  for  the  normal  distribution  with  mean  1  and  variance 
of  .5  is  e(-*25)+2(-*25)2  (.5),  Hence,  the  first  design 

will  have  a  present  worth  of 

-•25+i(-.25)2(. 5) 

$100,000  +  $15,500  1 _ _  =  $158,685.79 

1_e-,25+z(-.25)2(. 5) 
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and  the  second  design  will  have  a  present  worth  of 

-.25+i(-.25)2(.5) 

$120,000  +  $10,000  Z _  =  $157,861.80  . 

1_e-.25+i(-.25)2(.5) 

Considering  the  funding  as  taking  place  at  random 
points,  the  advantage  has  changed  to  the  second  design. 

The  wrong  alternative  has  been  selected  by  using  deter¬ 
ministic  data.  In  addition,  if  one  assumes  that  the 
expected  true  cost  of  the  first  design  is  $158,685.79 , 
using  the  estimated  cost  of  $154,572.58  yields  an  error  of 
$4113.21.  Obviously,  this  data  was  selected  to  make  the 
point  that  both  incorrect  selections  and  incorrect  estimates 
can  be  made  when  data  is  used  in  a  deterministic  manner 
even  though  the  data  is  known  to  be  random  variables , 

The  variance  for  the  present  worth  is  found  by  using 
only  the  variances  of  the  ^sd-^ss  model  element. 

V(PW)  =  c2  [Mt(-2r)-(Mt(-r) )2] 

[l-Mt(-2r)  ]  [l-Mt(-r)j  2 
where  C  =  $10,000 

M+(-2r)  =  e-2(-25J+2C-.25)2C-5) 

X 

=  .6456 

=  e-(-25)+i(-.25)2(.5) 

=  .7911 

V(PW)  =  $128,416,795.60 
0"pw  “  $11,332.11 


It  should  be  noted  that  in  finding  the  expected  value 
and  variance  of  the  cash  flow  a  specific  distribution  of 
cash  flow  has  been  assumed.  As  a  result,  the  distribution 
of  the  present  value  can  be  assumed  to  be  normally  dis¬ 
tributed.  Since  the  variance  of  the  present  value  is  known, 
more  information  is  known  than  by  the  first  method. 

To  illustrate  this,  let  us  assume  that  the  analyst 
used  the  data  in  a  deterministic  manner  and  selected  the 
first  design.  Since  the  true  expected  present  worth  of  the 
first  design  is  $158,685.79,  one  question  might  be,  what 
is  the  probability  that  the  actual  cost  will  be  less  than 
or  equal  to  the  present  worth  of  $154,572.58  calculated  by 
the  analyst?  Assuming  present  worth  is  normally  distributed 
with  mean  $158,685*79  and  variance  of  $128,416,795*80, 

Z  -  X-//  _  154,572.58  -  158,685.79  _ 

a”  ~ -  ~  -•36* 

C>  11,332.11 

Hence,  from  normal  tables  (Miller  and  Freund,  1965*398) 

(X- •359'  As  a  result,  there  is  only  a  35*9  percent  chance 
that  the  first  design  will  have  an  actual  cost  less  than  or 
equal  to  the  cost  estimated. 

In  many  applications  of  this  research,  multiple  cash 
flows  will  be  seen  in  the  models.  As  a  result,  it  is  possible 
that  no  distributional  assumptions  concerning  the  random 
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variable  of  present  worth  cam  be  made.  In  such  cases  it  is 
still  possible  for  certain  probability  statements  to  be 
made  via  the  Ghebyshev  inequality.  Using  the  Chebyshev 
inequality  for  this  example,  the  probability  of  realizing 
an  actual  cost  within  two  standard  deviations  from  the 
mean  is 

P  (  jPW-E(PW)|  2  (f  )  =~  l-i_  =  75#. 

k^ 

Another  way  of  stating  this  is  that  the  probability  of 
being  within  20,($22, 664. 22)  of  the  expected  present  worth 
is  greater  than  75#. 

The  strength  of  Chebyshev' s  theorem  is  that  it  only  needs 
the  mean  and  standard  deviation  of  the  distribution.  It  is 
also  its  greatest  weakness  in  that  it  only  provides  an 
upper  bound  for  the  probability. 

Other  tests  could  be  made  if  we  made  distributional 
assumptions?  however,  the  actual  distribution  is  not  avail¬ 
able  to  test  our  assumptions.  The  suggestion  is  to  then 
generate  the  distribution  via  computer  simulation,  validate 
the  expected  present  worth  and  variance  of  present  worth 
via  the  analytical  tools  given  in  this  research,  and  then 
make  probability  statements  after  analysis  of  the  distri¬ 
bution  of  present  worth. 

The  next  example  is  a  replacement  problem  first 
introduced  by  Grant  ( 19 38 *  208 ) .  The  problem  is  revised  to 
introduce  random  variables  into  the  problem. 


Replacement  Economy 

In  replacement  economy,  an  in-service  piece  of  equip¬ 
ment  is  considered  for  either  replacement  or  improvement. 

The  factors  which  govern  the  need  for  thi.<  nalysis  is 
obsolescence,  inadequate  capacity,  deterioration,  improved 
equipment  availability,  etc. 

The  following  example  is  a  decision  process  in  extending 
the  life  of  an  item  versus  immediate  replacement. 

Example  7*2  A  wooden  telephone  pole  has  decayed  to 
the  point  that  the  pole  must  be  either  removed  or  "stubbed". 

A  stub  costs  $8.00  with  an  installation  cost  which  is 
exponentially  distributed  with  a  mean  of  $20.00.  The 
inspector  notes  that  the  upper  section  of  the  pole  should 
be  good  for  about  5  years.  He  estimates  that  this  estimate 
is  normally  distributed  with  a  mean  of  5  years  and  a 
standard  deviation  of  1  year.  A  new  pole  will  cost  $25*00 
with  an  installation  cost  which  is  exponentially  distributed 
with  a  mean  of  $40.00.  A  new  pole  has  a  normally  distri¬ 
buted  life  of  20  years  with  a  standard  deviation  of  4  years. 
The  company  has  an  existing  contract  to  buy  all  used  tele¬ 
phone  poles  at  the  rate  of  $5.00  per  pole.  The  minimum 
required  rate  of  return  is  23  percent. 
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Figure  7.2  illustrates  the  decision  process.  Several 


model  elements  are  used  in  the  analysis.  The  initial 
procurement  costs  of  the  stub  01  new  pole  is  a  constant 
and  is  simply  added  to  the  present  worth.  To  calculate 
present  worth  due  to  installation  costs  and  the  salvage 
value,  two  models  presented  in  this  research  are  required. 
Calculating  present  worth  due  to  installation  costs  will 
require  the  use  of  the  CdgTcd  model  element.  The  CddTd 
model  element  is  required  to  calculate  the  present  worth 
due  to  the  salvage  value . 

The  expected  present  worth  of  the  stubbed  pole  is: 

E(PWS)  =  8  +  5  +  20(1)  -  5(e"23(  5)+l(--23)  (D2) 

=  $31-37  . 

and  the  expected  present  worth  of  the  new  pole  is: 

E(PWn)  =  25  +  40(1)  -  5(e_*23(20)+i(-,23)2(M2 
=  $64-92  . 

It  will  apparently  pay  to  stub  the  old  pole.  The 
difference  in  expected  value  was  really  determined  by  the 
high  cost  of  installing  a  new  pole.  Using  the  means  of 
this  data  in  a  totally  deterministic  manner  would  result 
in  the  same  conclusion  with  very  close  figures  to  those 
calculated  when  using  the  data  as  random  variables.  The 
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main  advantage  in  using  random  variables,  in  this  case, 
is  that  one  can  now  calculate  a  variance  and  standard 
deviation  for  the  estimated  present  worth  of  the  stubbed 
pole  i 

V(PWS)  =  (20)2  +  j2  [e-2(.23)5+i(-2(.23»2(l)2] 

.  [e-<*23)5+i(-.23)2(D2]2 

=  400.14 
PWg  =  $20.00 

Having  the  standard  deviation  does  give  one  a  great 
deal  more  information.  If  one  uses  deterministic  data, 
the  present  worth  of  the  stubbed  pole  is  $31.42.  In  this 
case,  the  advantage  in  using  this  research  is  in  being  able 
to  obtain  variance  and  standard  deviation  data. 

The  next  example  is  a  lease  or  own  problem  taken  from 
Taylor  (1964i21Q).  The  problem  is  revised  to  introduce 
random  variables. 

Lease  or  Buy 

To  own  or  lease  equipment  is  a  common  problem  in 
economic  analysis.  The  obvious  advantages  of  leasing  is 
in  avoiding  many  of  the  costs  of  ownership  such  as  obsole¬ 
scence,  repair  maintenance,  and  replacement.  The  following 
example  is  applied  to  machinery;  however,  the  concepts  are 
applicable  to  many  investments  in  capital  equipment. 


Example  7 . 3  A  2-year-old  truck  has  a  net  realizable 
value  of  $3,000  and  is  expected  to  have  a  salvage  value  of 
$900  after  its  remaining  3-year  life.  The  actual  salvage 
value  is  normally  distributed  N(900,450).  Its  operating 
disbursements  for  taxes,  insurance,  and  registration  are 
$160  a  year.  Annual  inspection,  maintenance,  and  repair 
costs  are  estimated  to  be  composed  of  two  costs.  The  first 
cost  is  for  preventive  maintenance  and  is  estimated  to  be 
normally  distributed  with  a  mean  of  $100  and  variance  of 
$900.  The  second  cost  is  a  ramp  function  which  is  expected 
to  increase  at  a  rate  of  $50  per  year.  The  $50  rate  is 
exponentially  distributed  with  a  mean  of  $50.  This  second 
cost  is  primarily  due  to  breakdowns  with  the  seriousness 
of  each  breakdown  increasing  with  time.  The  time  between 
these  breakdowns  is  poisson  distributed  with  a  mean  of  one 
year.  An  equivalent  truck  can  be  leased  for  20  cents  a 
mile  plus  $15  a  day  for  every  day  that  the  customer  keeps 
it,  whether  it  is  driven  or  not.  The  annual  utilization 
cost,  based  upon  past  records,  is  expected  to  be  normally 
distributed  with  a  mean  of  $1,050  and  a  standard  deviation 
of  $500.  Minimum  required  rate  of  return  is  \yfo. 

Figure  7.3  illustrates  the  network  for  owning  the 
truck  and  the  network  for  leasing  a  truck.  At  node  2,  in 
both  networks,  the  variable  N  is  increased  by  the  constant  1 


Owned  Truck 


Leased  Truck 


No 


Figure  7.3 

Owned  vs.  Leased  Truck 
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to  designate  the  period  of  "current"  analysis.  The  input 
of  ”R(E(50))"  at  node  5  designates  a  ramp  function  with  a 
cost  parameter  which  is  exponentially  distributed  with  a 
mean  of  $50.  The  method  of  branching  is  tested  after  node  3 
and  5  in  the  "owned  truck"  illustration  and  after  node  2  in  the 
"leased  truck"  illustration.  This  test  is  on  the  "current" 
time  for  the  alternative.  The  actual  network  arrangement 
is  discussed  in  more  detail  under  Suggestions  for  Further 
Research  in  the  next  chapter. 

The  analysis  for  expected  present  worth  for  the  owned 
truck  is: 


E(PW0)  =  3,000  +(100+160) 


Mti(-r)-Ki(-r)]4 


+  50 

■  l+Mt2(-r)3  [(Mt2(-r))3-4]' 

[l-Mt2(-r)]2  J 

- 

where  M^_1(-r)  = 
Mt2(-r)  = 

Mt2(-r)  = 
Mt3(-r)  = 


900  [M  (-r)]  , 

e-15 

e(e",15-l) 

_e(e",15-1-.15) 

e-45  . 


15? 


The  resulting  calculation  yields  a  present  worth  of 
$3*195*90.  The  analysis  for  the  expected  present  worth  for 
the  leased  truck  is 


E(PWl)  =  1050 

*Mt(-r)  -  [Mt(-r)]4' 

l-Mt(-r) 

where  M^(-r)  =  e-*  0  . 

The  resulting  calculation  yields  a  present  worth  of  $2,351.11* 
In  this  case,  present  worth  is  dealing  with  costs  and  the 
leased  truck  is  the  cheaper  alternative. 

These  examples  have  shown  the  usefulness  of  this  research 
in  extending  traditional  or  classical  engineering  economic 
problems.  These  are  two  primary  benefits  that  are  seen  by 
the  author  in  this  effort.  First,  it  is  rare  that  one  can 
forecast  costs  or  profits  with  certainty.  This  research 
enables  one  to  use  a  more  realistic  approach  to  using  the 
known  data  in  a  problem.  The  second  benefit  is  to  the 
student  in  that  there  should  be  an  increased  motivation  in 
dealing  with  analytical  tools  which  approach  real-world 
problems  with  fewer  assumptions  and/or  simplifications. 

The  next  area  of  discussion  for  the  use  of  this  re¬ 
search  material  is  for  the  validation  and  verification  of 


simulation  models 


Validation  and  Verification  of 
Simulation  Models 
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One  of  the  drawbacks  in  attempting  to  solve  real-world 
problems  is  the  complexity  of  the  available  mathematical 
tools.  Another  drawback  is  that  appropriate  mathematical 
tools  may  not  be  available.  An  approach  which  can  be  used 
to  address  both  of  these  problems  is  simulation.  Also, 
simulation  of  an  alternative  can  produce  results  beyond 
the  analytical  calculations  of  mean  and  variance.  As 
pointed  out  earilier,  using  Chebyshev's  inequality  to  make 
probability  statements  can  at  best  produce  weak  probability 
statements.  Making  other  assumptions  such  as  normality 
will  produce  errors  if  the  present  worth  distribution  is 
abnormally  skewed  or  peaked.  In  using  simulation,  the 
actual  distribution  of  present  worth  is  available  to  make 
probability  statements. 

To  have  any  faith  in  the  output  results  of  a  simulation 
program,  one  must  be  able  to  validate  and  verify  the  model. 
The  analytical  tools  presented  in  this  research  should  aid 
the  model  builder  in  both  of  these  efforts  by  allowing 
validation  and  verification  of  expected  value  and  variance 
data. 

The  next  area  where  this  research  could  be  put  to 
practical  use  is  in  the  development  of  a  simulation  language 
for  engineering  economic  problems. 
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Development  of  a  Simulation  Language 

Simulation  is  a  relatively  new  technique  in  solving 
capital  investment  problems.  Hess  and  Quigley  (1963)  were 
among  the  first  to  use  Monte  Carlo  simulation  techniques 
for  the  construction  of  output  distributions  for  present 
worth.  Other  authors  which  have  used  simulation  include 
Hertz  (1964),  Bussey  and  Stevens  (1971),  and  Whitehouse 
(1974). 

When  dealing  with  series  type  cash  flows,  there  are 
two  basic  approaches  which  are  used  in  the  literature.  One 
approach  is  rather  brute  force  in  repeating  a  node  for  the 
number  of  times  required  for  the  series.  An  alternate 
method  would  be  to  "loop"  back  to  the  node.  Both  of  these 
methods  are  time  consuming  as  far  as  computer  time  is  con¬ 
cerned.  An  alternate  approach  using  this  research  is  to 
use  a  single  node  to  calculate  the  present  worth  of  a  series 
of  cash  flows.  As  an  example,  let  us  assume  that  both  the 
cash  flows  and  the  time  between  cash  flows  are  random 
variables.  For  a  specific  run  of  the  simulation  program, 
which  would  generate  one  value  for  present  worth,  two 
random  number  generators  take  a  sample  from  the  cash  flow 
and  time  distributions.  The  present  worth  for  this  specific 
realization  would  then  be* 

n 

I  e-rt  . 
t=l 


PW  =  c 


The  variable  n  represents  the  total  number  of  cash  flows 
allowed,  and  it  could  be  determined  by  dividing  the  economic 
life  by  the  time  parameter.  Si'  a  this  would  generate  a 
real  number  in  most  cases,  a  simple  procedure  would  be  to 
truncate  n  to  an  integer.  For  example,  if  the  time  parameter 
was  .3  years  and  the  economic  life  was  10  years,  n  could  be 
33.  Repeated  runs  through  this  node  would  produce  the  same 
expected  value  of  present  worth  and  variance  of  present 
worth  as  calculated  using  formulas  from  this  research. 

Using  this  type  of  a  node  will  allow  the  analyst  to  validate 
segments  of  the  model.  Since  the  actual  distribution  of 
present  worth  would  not  be  available,  the  actual  use  of 
these  nodes  would  be  in  building  a  model.  Once  the  model 
is  ''correct",  these  nodes  would  be  replaced  with  appropriate 
nodes  to  collect  data. 

In  this  chapter  the  main  objective  was  to  examine 
applications  for  this  research.  Two  basic  areas  were 
discussed,  one  area  being  the  direct  use  of  this  research  to 
analyze  problems  and  the  other  area  being  the  indirect  use 
of  this  research  in  additional  research  areas.  The  next 
chapter  is  a  summary  of  this  research  and  includes  suggestions 
for  further  research. 


Chapter  8 


SUMMARY,  SUGGESTIONS  FOR  FURTHER 
RESEARCH,  AND  CONCLUSION 

The  objective  of  this  research  was  to  develop  analytical 
tools  to  solve  specific  classes  of  capital  investment 
problems  in  engineering  economics  involving  mutually 
exclusive  alternatives.  A  total  of  16  basic  models  were 
developed  for  impulsive  and  series  cash  flow  model  elements 
under  the  assumption  of  deterministic,  constant  interest 
rates.  A  review  of  Table  3*7  indicates  a  great  deal 
additional  research  is  required.  This  chapter  summarizes 
the  research  and  makes  suggestions  for  further  research. 

Summary 

The  first  phase  of  this  research  was  to  develop  a 
detailed  classification  for  capital  investment  models.  A 
spanning  set  of  models  found  in  the  literature  were  reviewed 
to  determine  types  of  models  and  required  parameters. 

Next,  a  taxonomic  structure  for  capital  investment  problems 
was  developed.  This  area  of  the  research  extended  Fleischer 
and  Ward's  work  (1977)  to  include  a  distinction  between 
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models  involving  one-time  cash  flows  (impulsive)  and  series 
cash  flows.  A  further  extension  in  the  taxonomic  structure 


was  made  by  looking  at  independent  and  dependent  models. 

A  short  form  classification  was  also  developed  in  Chapter  3* 

A  total  of  sixteen  basic  models  were  developed  for  impulsive 
and  series  cash  flow  model  elements  under  the  assumption  of 
deterministic,  constant  interest  rates  (see  Tables  4.4-4.11). 
Dependent  model  elements  were  extended  to  allow  for  a  finite 
number  of  cash  flows  (see  Table  4.12).  In  addition,  a 
general  model  for  analysis  was  developed  in  Chapter  3  and 
expanded  in  Chapter  4.  Chapter  5  is  an  indepth  discussion 
of  the  contribution  of  this  research.  Chapter  6  introduced 
a  network  logic  to  modeling  economic  problems  and  Chapter  7 
discussed  application  of  the  research. 

This  research  has  highlighted  various  areas  where 
additional  research  is  required.  The  next  section  of  this 
chapter  is  a  review  of  these  suggestions  and  the  introduction 
of  other  research  areas  not  covered  in  this  research. 

Suggestions  for  Further  Research 

Although  a  great  deal  of  work  has  been  accomplished 
by  many  authors  in  the  areas  of  simulation,  risk  analysis, 
utility  theory,  etc.,  a  large  void  still  exists  for  many 
of  the  basic  models  which  are  required  in  analysis  of 
economic  problems.  Further  research  is  required  in  the 
development  of  these  models  and  the  use  of  these  models. 


Non-constant  Interest  Rate  Models 


As  pointed  out  in  this  research,  with  the  exception  of 
Reisman  and  Rao  ( 1973 ) »  stochastic  interest  rate  models 
have  not  been  developed.  The  argument  for  considering  the 
interest  rate  as  a  constant  is  founded  upon  the  concept  of 
a  minimum  attractive  rate-of-return.  Currently,  interest 
rates  are  changing  rapidly.  There  is  reason  to  believe 
that  management  may  not  be  able  to  either  determine  a 
minimum  attractive  rate-of-return  or  must  use  a  very  high 
rate  to  cover  risk  in  the  analysis  process. 

It  would  seem  that  an  area  for  additional  research 
would  be  to  develop  models  for  stepwise  constant  and  variable 
interest  rates  with  stochastic  parameters.  As  models 
become  more  and  more  complex,  simulation  becomes  an  impor¬ 
tant  tool;  however,  analytical  tools  are  still  required  to 
validate  and  verify  the  simulation  output. 

Alternate  Decision  Criteria 

The  decision  as  to  the  best  alternative  has  been 
assumed  to  be  based  upon  the  alternative  with  the  minimum 
expected  present  worth.  Other  decision  criteria  need  to 
be  researched.  Annual  cost  analysis  is  commonly  found  in 
economic  analysis.  It  would  not  be  too  difficult  to 


extend  this  research  to  introduce  this  decision  criterion. 


ft 
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Other  criteria  such  as  the  aspiration  level  principle  and 
the  most  probable  future  principle  should  also  be  explored 
for  use  in  the  analysis  of  data  generated  by  the  use  of  the 
analytical  tools  developed  in  this  research. 

Introduction  of  Additional  Parameters 

The  basic  parameters  used  in  this  research  were  first 
cost  (P),  revenue  (R),  expense  (E),  and  salvage  value  (S). 
Other  parameters  such  as  costs  due  to  taxes  and  profits 
due  to  depreciation  could  also  be  introduced.  One  method 
to  introduce  these  parameters  is  to  look  at  expenses  and 
costs  which  appear  in  a  given  year  and  calculate  an  expected 
net  cash  flow  and  variance  of  net  cash  flow  for  each  year. 
These  data  calculations  could  then  be  used  to  calculate 
the  expected  present  worth  for  the  alternatives.  The  data 
calculations  would  have  to  be  handled  as  n  impulsive  cash 
flows.  The  suggested  formulas  for  calculating  expecxed 

net  cash  flow  and  variance  for  net  cash  flow  are: 

j  k 

ECNCFi)  =  (l-00(  £  R.-  £  Ek)  +  OCCD) 

i=l  J  i=i 

and 

j  k 

V(NGF1)  =  (1-002  (  I,  V(Rj)  -  V(Ek))  +  «2V(D) 


Conclusion 


This  research  has  expanded  the  analytical  tools  avail¬ 
able  to  accomplish  the  analysis  of  economic  models  involving 
mutually  exclusive  alternatives  with  impulsive  and  series 
cash  flows  under  the  assumption  of  deterministic,  constant 
interest  rates.  The  use  of  these  tools  allows  the  analyst 
to  use  available  data  in  a  more  realistic  mode  and  calculate 
a  more  accurate  and,  in  some  cases,  more  correct  solution. 
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APPENDIX  I 

:d  mathematical  proofs 
and  derivations 


(1)  Proof  that  E(e"rt)  -  e“r 

Given*  T  is  a  random  variable 
To  prove,  E(e"r^ )  —  e“r  E<T> 

yjn 

Solution,  The  Taylor  series  expansion  for  e“r  about  the 
point  T  =  is, 

e-rT  =  e-rMt  -  r<T-/Zt )e“r^t  +  f2(T~ +  .. 

*  2i 

Taking  the  expected  value  of  both  sides  of  the  previous 
equation  leads  to, 

E(e“rT)  =  e-r^t  +  E(T-Mt)re“r//t  +  E(T“^t)2  _  + 

2» 

Recognizing  that 

E(T-Mt)  *  0 


e(t -Mt)2  =  CTt2  , 


2^2 


,-**)  =.  1  +  — 2-i  +  ... 

9 


Since  all  terms  in  the  expansion  are  greater  than  or 
equal  to  0,  the  right  side  of  this  equation  is  greater  than 


Since 


-r  E(T )  _  -r(l  t 
6  *“  “ 

E(e-rT)  :  e-rE(T) 


9 
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(2)  Derivation  of  the  variance  of  present  worth  for  C^dT^g 
Independent  Model 

Given i  PW  =  C  e“rT 

T  is  an  independently  distributed  random  variable 
C  is  a  constant 
V(PW)  =  E(PW2)-[e(PW)]  2 

To  find i  V(PW) 

Solutions 

(PW)2  =  C2(e“rT)2  =  C2  e_2rT 
E(PW2)  =  (C2) •E(e_2rT) 

E ( PW )  =  (C)*E(e-rT 

(E(PW ) ) 2  =  C2*(E(e-rT))2  . 

Hences 

V(PW)  =  C2  E(e_2rT)-C2(E(e"rT))2 

=  G2  [E(e“2rt)-(E(e"rt))2]  . 

Since 

E(e"rt)  =  M^(-r)  , 

(E(e"rt))2  =  (Mt(-r))2  , 

And 

co 

E(e-rT)2  =  Z  (e-rt)2  P(t) 
t=0 

oo 

=  I  <e-2rt)  P(t) 
t=0 

=  Mt(-2r)  , 


1 


A. 


A. 


Jm 


1 79 


Then 


V(PW)  =  C2  [  Mt(-2r)-(Mt(-r))2] 
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(3)  Derivation  for  the  mean  and  variance  for  the  present 
worth  for  a  CsdTss  model  element. 


E(PW)  derivation  ^ 

OO  y 

V  r.  ^  ti 

Given*  PW  =  C  L.  e"ri=l 

0=1 

t^  are  independent,  identically  distributed 
random  variables 
C  is  a  constant 
To  find*  E(PW) 

Solution* 

PW  =  c  L  e'ri=l 
0=1 

oo  0, 

E(PW)  =  E(cX  e“ri£ltl)  . 

0=1 

Since  C  is  independent  of  t^  , 

oo 


E  ( PW )  =  E  ( C )  •  E  ^  e 

0=1 

Since  C  is  a  constant, 


•E(  E  e-i^S 


OO 


0 

r"1  I 

E(PW)  =  C *E(  L  e"ri=l  ). 
0=1 


By  expanding, 

E ( PW )  =  C  [E(e"rtl)+E(e"r^tl+t2))+E(e“r^tl+t2+t3))+. . .] 


«-»»•  -T  ^E[?^SSSSOLUTlrNHT«^^s”EoSS«vrMlK  ECONOMIC  ,N»L«.S_ETC  <U> 
DEC  81  J  H  ESTES  HU 

UNCLASSIFIED  AFIT-CJ-81-66D 
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For  all  i  and  k, 

E(e“rti)  ■  E(e“rtlc)  «  E(e_rT)  . 

Setting  Efe"**1)  «  X, 

oo 

E(PW)  «  C  I  Xn  . 
n«l 

Since  le”1^  I**  1  for  all  T-'O,  X  1,  and  this  geometric  series 
converges  to 


Hence , 


Using 


E(PW)  »  C 


E(e_rT) 

l-E(e"rT) 


E(e~rT)  =  Mt(-r)  , 


E(PW)  «  C  • 

l-Mt(-r)  . 

To  findi  V(PW) 

By  definition  of  variance* 

V(Ptf)  a  E(Ptf2)- [E(PW)]  2  . 
By  squaring  PWt 


.  C2(.-rt*  ♦  .-r(t‘+t2>  ♦  .*r(t»*W  + 

•<.-«!  ♦  ,-r(tj«t2)  +  +  ...) 

«  C2(.*r2tl  ♦  .-r(2tl+t2)  ♦  ♦  ...) 

♦  c2(.-r<*t1»t2)  ♦  .-r(»i.2t2>  «.  #-r(*1*2t2*tJ)  +  , 

♦  c2(e'r(2tl**’tZ4,t3>  ♦  #-r(2tx4-2t2*t3>  «.  #-r(2ti«>2t2*ttj) 

♦  ♦  ...)  . 

Taking  tha  expected  value  of  both  aides  and  letting 
X  »  E(e*’2rt)  and  Y  ■  B(*“rt)  leads  to* 

B(PW2)  -  C2^  [  X  ♦  XY  ♦  XY2  ♦  XY3  ♦XY4  ♦  ...  ] 

♦  [xY  ♦  X2  ♦  X^  ♦  X^2  ♦  xV  ♦  x^  ♦  ...] 

♦  [xy2  ♦  X*Y  ♦  x3  ♦  xh  ♦  x3y2  ♦  x5y3  ♦  xV  ♦  . . .] 

♦  [xy3  ♦  x2*2  ♦  x3y  ♦  x4  ♦  x*hr  +x*hr2  ♦  x*y3  ♦  . . .  ] 

♦[•••]]• 


oo  oo  oo 

I  i‘  I  tU  I  Is  I  rJ 

k»l  J«0  k»l  J*1 

C2[  x  1  ♦  X  Y  1 

[  T^X  '  T-7  TT  *  I^Y  J 

p2  f  X*XY  l.„2  X(  1»Y )  . 

I  (l-X)(l-Y)J  (1-X)(1-Y) 
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Mow,  substituting  X  •  ^  *  »  B(s  *^1 

«*» .  c*  « (uet<:^ . 

noting  E(«"2rt)  •  HtC-ar)  and  Ets'^)  •  *t<“rK 

BO*2)  •  C2  -1— - 7 — -  • 

Ci-H^-EriHl^-r)) 

Sines  (fro*  p.  181) 

*A-r) 


ECPV)  •  C 


1-^t-r) 


_  f  *,(-r)  I2 

<B<**))2  •  C2  — — — •  . 

(.l-i^C-ry 

Hones. 

s(p»)  •  BCP^MtCP*))2 

.2  *s<-2r>t>#*v»r)) 

( l-l^t-ErW  l-&t(-r)) 

.  ^  j^j* . 

ll-K^-r)' 

which  can  bo  reduced  to 


•  C 


t<PU) 


•  C2f  Xt<-*rMHt(-r))* 

['l-a^-ar)]  [!-**(. r)p. 


i 


1S4 

(%)  OirintUn  for  tho  vtrUnet  for  tiw  prooont  norU  for  » 
CMf|f  ««M  oiaa  out. 

Cl  von  i  W  «  C|#'rtl  ♦  >  ... 

t|  art  independent .  i district** 

rondo*  nriiblH 

«l  oro  lAdfpmdMi.  disirstwvod 

rwdM  vtrUUn. 
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To  flndi  ft rr) 

MttUoni 

fey  definition 
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ly  oqooriAft  f»« 

♦  6  ,C jO-r^  #  #jVrtn,,}1,' 

♦  »  ...  *  CjOj* 

o  *  «>VHnl#n«*n^  •  ... 

Taking  uw  ecpeeted  vale*  or  toll  tldoo  and  lotting 

X  *  Ito-art)  and  T  *  •to*rtl#  <oada  to. 

tire*)  *  stc/a^stcjC^JCT^etcjCjitr3*... 

♦  MojOjW^tta^BjU^retldj*)*5*... 
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.[?  (SI 

Ln*2 

,  c — 

‘  i2  *  4  l 


e(t«'2rt)  •  B(te'2rt>MEU'2rt))n'2 


1*1 

wfcicn  is  equal  toi 
'  •  ^  ( «2r ] 7 • 


^  2  <*  *  s 

*  [  i  i  :  1  '•  2[  1  1 

U-^l-2r)l  U-Kt(-r)i 

*  U>  *  f?  ♦  4  f  \{mTi  ], 

'  *  !-*.(-«•)-- 


where 


V  tft*  •  ;r t  •  ?r*.  .  •  2rt  n-2 

*  *  i.  [,[  ttli*  ntSUe  ))*(E(o  )) 


■*> 


2  he  r-oztio?-  ft>r  a  can  no*  fee  written  as< 


A  *  (JJtt# 
<»o 


■*"»*  •  ?.  i; 

n»2 


1  ,  -2rt .  n-2 
I  (B(e  )) 


V  kM  -Crt  n.g 

/,  t'l  i££«  Ji  can  be  solved  using  the  following 


tips 


Le«  0*  E4e'2fl)  . 

then  fjQ!  *  )  >  pn  *  _L  . 

i-p 


f’fpl  *  i.  nPCn“n  *  1 


n*  1 


(1-pl 


2  * 
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and 


oo 

*"<P>  =  I  n(n-l)p 


n-2 


n=2 
n\  n-2 


(1-P) 


3  ‘ 


Since  „ 

L  Op 

can  be  written  as 


nH; 

oo 


,  x  n“2 

n(n-l)p  , 


I  ("I  (E(e-2rt»n-2 

n«2  ' 

A  can  then  be  written  ass 

.  ,  -2rt..2 

(E(te  )T  *  ' 


(1-E(e“2rt))3  ‘ 


(l-E(e-2rt))3 

and  can  be  solved  using  moment  generating  functions  ass 

2 


A  = 


-  Mt  (-2r) 


2  J  (l-Mt(-2r))3  ’ 


2 

Combining  both  sets  of  terms  (t^  and  1;^!:  * )  simply  leads  tos 


(PW2 

=  Mt"(-2r) 

1  +  2Mt(-r) 

l^r 

4(l-Mt-(2r))2 

l-Mt(-r). 

Mt' (-2r) *Mt(-r) 
(l-Mt(-2r))(l-Mt(-r); 


+  (Mt  (-2r))2 
(1-Mt(.2r)) 


3' 


1  +  2Mt(~r) 

l-Mt(-r)J 


- -j 
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(6)  Derivation  of  the  expected  present  worth  for  the  C__T__ 

So  SS 

dependent  ramp  model  element. 


Givens  PW  =  c1t1e~rtl  +  c^tj+tg)  e~r('bl+'t2^  +  ... 


-r(t1+t?+. . .tn) 


cn(t1+t2+...tn)  .  2 


•  •  •  • 


and  t^  are  independent  and  identically  distributed 
random  variables. 

To  finds  E(PW) 

Solutions 


Since  is  independent  of  tj_, 

E(PW)  =  E(o1),E(t1e‘rtl)  +  E(e2),B(t1e"r*tl+t2Vt2e'r(tltt2)) 

+  ...  +  E(on)'E(t1e-r(tl+^+---+tn)+t2e-r(ti+t2+...+tn) 

+  ...tne-r(tl+t2+"-«n,)+  ... 

Since  E(c.  )  =  E(c-)  =  fJL  for  all  i  and  j, 

E(PW>  ■=  fJLc  [E(t1e'rtl)+B(t1e'r<tl+i2,)+E(t2e"r(tl+t2))+... 

+  E( t^1 ‘l«2t-  •  •+t")  )+E(  t2e-r(tl+t2+-  •  •+tn) )+. . . 

♦  E(tne-r(tl+t2+”-+tn) )+...]  . 

Since  the  t^'s  are  independently  distributed  random  variableSi 
E(PW)  =  JUC‘  E(t1e“rtl)  [l  +  E(e'rt2)+E(e"r(t2+t3) )+...] 


LiA 


X 
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+  Mo-E(t2e'rt2)-E(e-rtl)[l  +  E<.'rt3)rt<.-rtV1*>>*...] 
+. . .+  McE(tne-rt")-E(e'rtl)  •E(e'rt2).  •  .zC^"  > 

[l  +  E(e"r<tn+1>)+E(e'r(tn+l+tn+2) )+•••] 

+  •  •  • 

However,  E(tie”rti)  =  E(t^e  rtj)  =  -  Mt 

and  E(e"rti)  =  E(e”rtj)  =  Mt(-r)  for  all  i  and  j. 


Hence i 

E(pw)  =  Mc  [-Mt'(-r 


)]  •  I  K<-r>]  "•  £  K<-r>) 

J  n=0  1  fc=0 


k 


=  Mc  (-Mt'(-r)) 

.  -  McMt,(;r)  . 
[l-Mt(-r)]2 


1 

l*Mt(-r) 


1 

l-Mt(-r) 


-  j - 
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(7)  Derivation  of  C88T88  Dependent  rasp,  finite  n  caah  flow 
model  element 

-rt,  .  ,  -rttj-ftj) 

PW  ■  Cjtje  1  ♦  c2(tj4-t2)e  1  2 

♦  e3(t1+t2+t3)e"r(tl+t2+t3>  ♦  ... 

♦  cn(t^t2+t3+. .  .+tn)e‘r(tl+t2+t3+*  * ' tn) 

-  cltle-rtl  ♦  c2t1e-r(tl+t2)  ♦  c^.'^VV 

.  -r(t1+t,+t*+...+t_)  ^  .  -r(t«+t2) 

+...+  c^t.e  l  z  3  n#  +  c~t~e  1  z 

n  i  z  z 

♦  c ,t,.-r'WV  c  t,.-r(tltVt3t"-+tn> 

3  z  n  z 

♦  c3t3.-r(tl+Vt3)  ♦  cnt3.-r(t»+t2+t3+",+t>>) 

+...♦  C  t  ,-«■< . . . nn) 

n  n 

and  t^  are  independent  and  identically  distributed  random 
variables. 

Solving  in  terms  of  expected  values i 

E(PW)  »  E(C)*E(te~rt)  +  2E(C),E(te"rt),E(e‘rt) 


+  3*E(C)*E(te’rt)*  [E(e”rt)]  2 
+  ...+  n  E(C)*E(te'rt)*  [E(e’rt)J  n  1 


n 


-rtjk-l 


E(C) *E(te"rt) •  Y,  k[E(e"rt)] 
k-1 


X  *  E(e-rt)  . 


Let 


Than 
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t 


k(E(# 


-rtnk-l  ,  x0  ♦  ♦  ytz  «* 


n-l 


.  2-  f  . 

**  & 


N0«» 


X  -  X**1 


?  x*.£-^ 

L  i  -  x 


k*l 


d 

—  S 
dx 


n 


(l-XKMn+nx")  ♦ 

U-X)2 

l.x  .  (tHl)Xn  ♦  Xtn*l)Xn  ♦  X  >  e"*1 
~  (1-X)2 

1  *  Xn(»n»l*RX»X»X) 

(l-X)* 

1  +  Xn(nX  »n«  1 ) 

(1-X)2 


n 


k»l 


f „  -rtlk-l 

1  ♦  [BC«"rt)]  ni.nSt*  r  )~n‘V 

[  E«  J  * 

<l-S(«'rt»2  ■* 

Sinc«*  S<w-”>  ■  -Kt'l-r)  • 

E(«'rt)  -  *tl-D  . 

and  E(C)  *  Mc 


R(PW) 


M  l^»1n  K^l 


1 


:?> 

{ *)  Derivation  of  Cg#?ss  de?*y,  fjr.jto  «  casn  f)t>* 

mo&s t  eleaon? 

rv  *  e.#*U*rUl  *  e*e 

*  V 


* i 1 j*1 

i4»»fH,.i,f*^******' 

•  •  C.  9  f  * 


'a 


1 9MtAg  •X?PC\*4  valua««  oAiio  W41  Mk»  fj'J  *n<S  T/S 

AT*  illoA*. rarfvOOJt  VATiAfclee, 

leads  10 

«<r»j  *  Ue  f 

*  iN  v  1 

Salting  Eie“***ri5}  *  Et 

Ze«.-‘**rV  ,  f.  **  *  *  t  *■>  .  *f  ud 

fc*l  **I  /*«  ,  J.*  ' 


Honcoi 


E(?V) 


•  i» 

ci  i  .  u«-U*rn 


). 


and  in  toraw  of  •o*ent  <?*nera;SA£  f ur*c  1  i on# • 


*ir*>  *  M, 


t  1  -  Kt{-(a*r)) 


n*l 


i 


i/ 1  i*  <5  o!  2  f;ui‘  •’■  <r&s.J. 

^C'  « C 


f  l  J*  flu^re  l  cl  on**-". 


•»*■*.  i  » t  "•  ■»  »  *  i  *.  *  *  1  ,  *  r  i  *.  v  •  i  .  1 

j*'#  *  ?  i  .  t  *  ;  c  *  »  -r  .  I  5 »«  ’  -  c  3 


1  •  %  •  *  4  *  .  •  *  ® 


T  -  - -■  t  *  r  T  • 


?ft*  i/f  o'v*. *  f  ♦;■■.«? ",  i'etf-  is  f  j*  ;  o’  *.*.c  sv«?i  ; ; t)‘.  «u  ■£  l.f.c 

4«ii  V/  f  «*•;-•  *.  ;  it--.  •  --*«  f  0)V.  iv4’  .  *1  *.  >JC  ;  c  *  »K  *.  I  #  i  c  j  <U  -t 

f*  ■  * 

i  «r*  *  i 

I  ^  .?•.  *  *  \  .  *»r  *,  vs. 

<£  l  #  ^  Li  L*  £  ■■ 9  **■  i  e  • 

: .  '  '  J 

4ft  $  iA-  j  {  £ «  *->*  f  *-.  ,  ;  »  .:4. 


*;?*:  *  14.  f  JT'  ‘  /-•  •"  . 


'  B  X  »*  t  •  I  •  / 


_  .  a«t 


.3:^'  -iiwW-t  i*  w*,i  t'iffc.  ;ju  t’l  i  ','*+A*»  o». 

l?.  I?*".  m»  **<*  «*** 

vft  8A*  fci  fv»ti#  &e44*ft4*.  }*,  jvij-  i* 

;»h.i  1 1#  i*-  *•«.  <• .  ;fw«M  *.’••■4  4«n  «)U«»W1 

■.«*mat*«l*^  »-i*  U-  }i<:.  -,-4*  «r4»#i*<  U<  *.4.* 

*4f  fr'tfe*  i*.  i>«<9  *#4  w»*  **)«£- 1*4  *^4*#'  t*,*  Aim**. 

ftfectsiA*  #*<*  *  <wwki #*  4  n*\u-+.  fragr**  i«  l*fc<  \*  «n«t4 

**l  *■**■■+  *.«.(«*  r<*4Vf  u-  ?«4»S*.  JiMiftt#.  ft*  r  «*$ 

Pi*  «f  fcei***«  ##£#"4*  *1  tfUMk  iut«  U4««ft;Ty 

Iip  ti*;  4*4  f*««i**«  M*  e«MU«|*ii.  M  «  £*e«*4  \ 

4«.  fit#  ■>. J?»’4  imi**  44 f  fm#  ijfe*  «4n#  3*4  *irtT» ntitr 

*f  l»7;  **  ft***?**  ATI*  #*»*  it***  U«i» 

«4SH»  N*u4#f  8A<*  4AT  fw*  ?«e*r-*3*Or  1* 

«  *i  ui*+e*  4+€T9+  l*  iMb««tri»t  iyp*ti**»i*r5*tf. 

8a*  i»?  *44  ;.M’,«rt  a#*#*.:  4*.  * *>***!> 

if  t p.#  v*c  #**4 ^404  *,*»  *if  fort#  *****  4*,  c*j j;^n,;*. 

''+**•  v-<4  : »;  ^i-*4w.  4*  *.-•  3r«f<#*4or  «f 

45  ?."!•  r.ft.i»,#4  £i.4i»«  Air  r«rt«  »*-  C«3»r*<«» 

AO  #4*  #•;•**!«.*  l*  rwioTv  U»  AT  4  &<>*■*  4t«1o  «*4*«r&i1#  It, 

t»"*l  84  r>4o«#  4  <»»TN  tf  s*e|«f  *f  *r  4# 

4  :##5*»ft.  1*  14*  siNiM  iui**  Air  fore#  «i4  u 
«  «f  :■*»  irt*  t|  •  j»v. ’.iw*rs**  *•*<»  3«ei#iv^  Alt** 

?i  t  iai  4#*e  *#*•■«.#  SotivljpJ.  Air  rorf# 

A*#4#  «it  «^ftr«  4T-4  A4i#r  4AW  #**ti  l**V*  Of  lt*4«s,Tt  i  A, 

*«**4#M  ft*  4*  tttd  «*#  riW  f*  *;  *#•*  #{#, 
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